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PREFACE 

When, man first invented the wheel, he did so without any theory. 

Several thousand years later, when man first went to the moon, he 
built the wheel for the Lunar Roving Vehicle (LRV) , again without 
adequate theory. If man had waited for exact theory, the LRV probably 
would not have been carried to the moon. 

So it seems that the existentialist is right in saying that 
existence precedes essence; that discovery proceeds from conceiving the 
physical phenomenon, to supplementing this phenomenon with explanation, 
usually expressed by mathematical abstractions. At least in wheel-soil 
interaction, wheels have preceded theory. Thus, a theoretical framework 
for a physical phenomenon which is known to exist can be useful, but it 
is not essential. 

This philosophical argument presents a fundamental restriction to the 
value of any theory. There are, however, other less fundamental but more 
immediate restrictions. 

In reviewing trafficability literature and in contemplating such 
phenomena during the past two years, the author (schooled primarily in soil, 
mechanics) has become impressed with the complexity of the general problem ot 
off-road mobility (see, for example, Bekker, 1969). This problem does not 
consist merely of a single, ideal wheel operating in an ideal environment, 
but of coupled wheels operating both in ruts made by preceding wheels 
and in virgin terrain. Further, the traverse over such a terrain is 
likely to include various soils, rocks, and bumps. The dynamic overall 
ride over such a terrain may, therefore, be only remotely related to 
theory developed for ideal wheel-soil conditions. 
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A general solution to a complex problem often results from research 
on separate aspects of the problem. This report presents a different and 
perhaps practically applicable theory on the very limited phenomenon of 
a rigid cylindrical wheel operating in a homogeneous soil. The develop- 
ments are based on relatively simple considerations of statics and dynamics. 
Fundamental observations render the problem determinate. This leads to 
solutions of the sinkage and the pull which are likely to be within 15% 
of the correct value. 

It is hoped that the theory presented in this report will be useful 
in evaluating and designing wheels for off-road mobility. Some new ideas 
are presented, and these need to be thoroughly checked and tested. It is 
hoped that concepts such as 1) the line of action of the resultant of 
radial stresses, 2) slip at a point, 3) the shear stress surface 
T * f(0,s), 4) the closed -form approximate relation between contact angles 
and sinkage, 5) the general graphical solution for pull, 6) soil inertia 
forces, 7) equivalent cohesion, and 8) the performance surface, will inspire 
further thinking and relevant research. 

The author wishes to express his appreciation to the National 
Aeronautics and Space Administration (NASA) for sponsoring this research „ 
Special thanks go to Professor James K. Mitchell for his continued support 
and patience throughout the investigation, and for reviewing the manuscript. 
The author also greatfully acknowledges Mr. Dieter J. Schuring for 
reviewing Chapter 7 and Ms. Madeline Travers for editing the report. 

Finally, in one way or another, this treatise reflects understanding and 
experience accumulated in off-road mobility literature, without which 
this study would have been impossible. 
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CHAPTER 1. INTRODUCTION 

A primary objective in off-road mobility is the design of a wheel 
which allows the most efficient transformation of mechanical energy of 
the engine to translational capability of the wheel. The most important 
single parameter in such a design is the pull, P, that the wheel can 
develop. This is defined as the pull that can be developed from the 
traction of a given wheel on a given soil for a certain input torque, T. 

In order to perform the desired design, it must be possible to predict 
this pull, either from theory or from experience. 

Previous approaches to wheel-soil interaction have been either 
quasi- theoretical (Bekker, 1956) or empirical (U. S. Army Waterways 
Experiment Station [WES], 1954). One major shortcoming of the empirical 
approaches, and to a certain extent of the quasi- theoretical approaches, 
is the necessity to test a full size wheel. Adequate theory would consider- 
ably reduce the need for such expensive testing, although it would not do 
away with testing entirely. 

A fairly substantial body of literature on wheel-soil interaction 
has accumulated over the years. An adequate review will not be attempted 
1) because an excellent review is available (Bekker, 1969), and 2) because 
most of the previous approaches have little direct bearing on the develop- 
ments presented in this report. 

Because the wheel-soil interaction problem is indeterminate, and the 

torque input and wheel variables are rather specific, many investigators 

* 

hoped that the secret to the solution could be found in the soil. Conse- 
quently, much recent research has dealt with methods of measuring and pre- 
dicting soil flow and pressure distributions at the wheel soil contact. 

7 

These investigations have resulted in a better understanding of the 
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wheel-soil interaction problem. They have also led to the conclusion that 
deformation processes in the soil are very complex, and that consequences 
of assumptions are difficult to evaluate. 

This report presents both rigorous (general) theory and approximate 
theory for wheel-soil interaction. The purpose of the general theory is 
primarily to form a basis for further developments. The approximate 
theory forms the basis for a practical solution to the problem. It will 
be shown that two fundamental observations render the problem 
determinate: 

1. The line of action of the resultant of radial stresses 
acting at the wheel soil interface approximately bisects 
the wheel-soil contact angle, 6 T , for all values of slip, s. 

2. A shear stress surface, T = f(0,s), can be hypothesized. 

The influence of soil inertia forces is also evaluated. A concept 
of equivalent cohesion is introduced which allows a convenient experi- 
mental comparison for both cohesive and frictional soils. This theory 
compares favorably with previous analyses and experimental data, and 
shows that soil inertia forces influencing the motion of a rolling wheel 
can be significant. It is assumed that the inertia forces can be super- 
imposed on forces resulting from a static analysis. 
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CHAPTER 2. WHEEL-SOIL INTERACTION 

As a wheel moves on a deformable soil, a complex interaction takes 
place at the wheel-soil contact and in the adjacent soil. Fig. 2-1 
illustrates a track and associated shear surfaces resulting from pulling 
a 24-inch (61-cm) diameter rigid spherical wheel on Yuma sand. 

Model studies of this interaction (Hovland and Mitchell, 1972) show 
that the most noticeable features are: 

1. Forward bending (movement) of initially vertical soil sections; 
some associated lateral movement. 

2. Volume change in soil (compression directly under the wheel and 
dilation to the sides and front) . 

3. Shear surfaces, and sliding primarily forward along these 
surfaces . 

To explain some of the details associated with general 
shear , the sequence of rolling and shearing is 
diagrammed in Fig. 2-2; (a) the soil before the wheel has rolled over 
it; (b) the shear surfaces and deformations that develop as the wheel 
rolls to the right, to position B; and (c) the appearance of the section 
after the wheel has rolled past. 

Observations of particular interest are: 

1. Shear surfaces do not appear to originate at the wheel surface 
(see Fig. 2-2b) , as is often considered the case with shear 
under a loaded footing. 

2. At the time a shear surface develops, its lower end is roughly 
parallel to the wheel surface. It diverges from a direction 

parallel to the wheel surface as it proceeds up and forward, 

> 

k wedge of soil is apparently pushed up and forward. 


3. 



Fig. 2-1. Track and shear surfaces resulting 
from pulling a 24-inch diameter, 
564-pound sphere on Yuma sand. 
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4. Movement along any one shear surface is restricted, as implied 

, by the relatively short but constant displacements. 

5. Shear surfaces are spaced at relatively constant intervals as 
can be seen from the spacing of the teeth marks (compressed 
remnants of once active shear surfaces) . 

It is possible to conclude that the shearing process is not 
continuous, but consists of separate small shear phenomena spaced at 
relatively constant intervals with respect to space and time. The 
spacing of the shear surfaces as well as the amount of movement along 
any one shear surface is undoubtedly a function of soil type. Hence, 
although the rolling of a wheel may appear continuous to the naked eye, 
it is possible that the soil acceleration acting on the wheel changes 
slightly from one shear surface to the next. 

The fact that shearing along any one shear surface appears to be 
restricted to a relatively small distance is perhaps associated with 
the change in direction of the shear surface with respect to the stress 
causing the movement. Most of the movement along a shear surface 
probably takes place when the surface is first formed. At that time, 
the shear surface is directed forward at an angle of 45° - (j>/2 to the 
direction of the major principal stress, As the wheel moves forward, 
the shear surface bends or turns toward a more vertical position. As 
the wheel passes, the shear surface is again bent down toward a final, 
more horizontal position (Fig. 2-2) . 

An instrumented spherical wheel, which was tested in Yuma sand, 
allowed further insight into wheel-soil interaction. The maximum radial 
wheel-soil contact pressure increased with increasing wheel load and then 
remained approximately constant, as showli in Fig. 2-3. This suggests 
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that the maximum radial pressure increases until the soil bearing 
capacity is reached. With further increase in wheel load, equilibrium 
is established by sinkage into the soil, with the load being distributed 
over a larger area. It is to be noted that, with the two lightest wheel 
loads, no distinct shear planes could be detected. Similar information 
can also be deduced from analysis by Vincent (1961) who noted that the 
rear wheel-soil contact angle for wheel loads too light to generate 
shear is relatively larger than the same angle for wheel loads heavy 
enough to generate shear. 

The number of shear surfaces that could be distinguished in the 
forward shear zone (soil in front of wheel) were counted, and the results 
are presented in Fig. 2-4 as a function of track depth. The point on the 
horizontal axis of Fig. 2-4 corresponds to the 1200-Newton wheel load in 
Fig. 2-3. 

The above discussion and data are primarily applicable to a spherical 
wheel rolling in sand. For such a soil, the interaction mechanism is 
dependent on the state of compaction. In a very loose state, the material 
would be compressible and the deformations would consist of both volume 
change and shearing distortion, but no general shear surfaces would develop. 
In a very dense state, general shear would be more significant, although 
volume change and shearing distortion would also take place. The usual 
case would involve all three phenomena; initial volume change accompanied 
by shearing distortion would be followed by general shear. 

The following steps appear to be involved as a wheel rolls over a 
sand surface: 

1. Compression occurs under the wheel. 

2. Initially vertical soil sections are bent forward. 

f 
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RESULTANT OF WHEEL LOAD AND TOWING FORCE IN 

NEWTONS 

Fig. 2-3. Maximum pressure vs. resultant of wheel load 
and towing force. 



Fig. 2-4. Number of failure planes within the forward 
shear zone as a function of track depth. 


/ 
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3. A horseshoe-shaped zone of dilation, which extends from the 
sides and around the front of the wheel, develops. This zone 
moves forward with the wheel. Some of the external evidence 
of this zone of dilation is the bow wave that forms in front 
of the wheel. 

4. When volume changes and shearing distortion can no longer account 
for all the soil that must be displaced, shear surfaces develop. 

The development of shear surfaces is probably influenced by the 
magnitude of shearing distortion and the magnitude and direction 
of the major principal stress. 

5. Sliding along the shear surfaces continues only for a short time 
and distances, as previously described. 

6. After a certain distance, a new shear surface develops and the 
shearing cycle repeats itself. 

7. As the wheel rolls forward, deformations assume a final position, 
and excess material from the bow wave is wasted to the sides to 
form the crests of the track. 

Other experimenters have investigated wheel-soil interaction in clay 
(Yong and Webb, 1969; Yong and Windisch, 1970). These studies indicate 
that deformations are similar in clay although the elastic recovery or 
rebound is greater. 

However, the relative importance of cohesion or friction depends on 
the nature of the phenomenon being investigated. A soil may be classified 
cohesionless by ordinary soil mechanics criteria, and yet for very small 
loaded areas cohesion may account for most of the resistance. Consider, 
for example, spheres rolling in a soil typical of lunar soil, with cohesion, 
c, of 20 psf (1 kN/m 2 ) and a friction angle, <p f of 37°. For analysis 
based on bearing capacity theory (Fig. 2-fp , cohesion provides most of 
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Fig. 2-5. Relative contribution to total soil resistance 
from cohesion, surcharge, friction and density, 
and soil inertia as a function of sphere diameter. 
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the resistance for a 10-cm in diameter sphere, while friction provides 
most of the resistance for a 2-m in diameter sphere. Fig. 2-5 also shows 
the relative contribution to total soil resistance (bearing capacity) 
from the various terms in the bearing capacity equation (Hovland, 1970), 
and it shows how these are influenced by the gravity field. 
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CHAPTER 3. GENERAL WHEEL-SOIL INTERACTION THEORY 

The developments in this chapter form a basis for theory in subse- 
quent chapters. It is hoped that these developments will also be of value 
to other wheel-soil interaction investigations which depend on a knowledge 
of the distribution of normal and shear stresses along the wheel-soil 
interface. 

DYNAMIC EQUILIBRIUM 

A free-body diagram of a wheel and all the forces and pressures 
acting on it are shown in Fig. 3-1. The symbols are defined below: 
a = linear acceleration of wheel 

u = angular acceleration of wheel 

m = mass of wheel, W/g 

W = wheel load or weight 

g = acceleration of gravity 

v = linear velocity of wheel 

w = angular velocity of wheel 

I = mass moment of inertia of wheel 
T = input torque 

P = output pull, pull developed by the input torque and thfe 
wheel-soil traction 
dN = differential normal force 
dF = differential shear force 

dR = differential resultant 

F a - algebraic, sum of shear forces along the wheel-soil contact 

R & = algebraic sum of resultants dR 

R - vectorial sum of resultants $R 

r = wheel radius 
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r^g = moment arm of any dR 
r^' = moment arm of R 

x,y,z = coordinates 

0 = any central angle measured counterclockwise 

8^ = angle between the vertical and front soil contact 

0£ = angle between the vertical and rear soil contact. 

For a rigid wheel rolling on an unyielding surface, the forces 

acting on the wheel can be represented graphically as shown in 
Fig. 3-2a. Analogously, for a wheel rolling on a deformable surface, 
the forces can be represented as shown in Fig. 3-2b. It is, therefore, 
valid to express the developed pull by 

P = ^ R 2 - W 2 (3-1) 


For the general case represented by Fig. 3-1, 



(3-2) 


The unknown in equation (3-2) is the vectorial resultant of all forces 

acting at the wheel-soil contact, R, which is a function of the shear 

and normal stress distributions. 

At any point on the contact surface, the normal stress may be - 

defined as a„ and the shear stress as 
0z 


In = c Q + a„ tanS Q 
0z a0s 0z 0s 


(3-3) 


If we further include the possibility of a toroidal, or other than a 

cylindrical wheel, the radius may also be a function of z. (Such cases 

have usually been studied using the average wheel radius called the 

effective radius, r .) Then the total shear force along the wheel-soil 

e / 

contact is for the general case 
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a) b) 

unyielding surface yielding surface 



Fig. 3-2. Force polygons of forces acting on a wheel. 




(3-4) 


F a = JJ r zKes + a 0z tan(S 0s )dedz 

In the above equations, the subscripts 8, s, and z indicate that the 

parameter having the subscripts is a function of the contact angle, 9, 

slip, s, and the dimension perpendicular to the plane of the paper, z. 

Also, in these equations, c & = wheel-to-soil cohesion or adhesion and 

6 = the wheel-to-soil friction angle. The mobilization of c and 6 

a 

will be considered later. They are in general much lower than the soil- 
to-soil cohesion, c, and the soil-to-soil angle of shearing resistance, cj). 

Taking moments about the center of the wheel (Fig. 3-1) and noting 
that all the normal forces go through the center, 

Rrl - F r = T - Iu (3-5) 

I cL 

Thus we note that the quantity Rr^’ is determinable if either the total 

shear force, F , or the input torque, T, and the dynamics of the wheel 
d 

are known. In this case. 


R = -V (T - Iu) (3-6) 

r f 

If the input torque and wheel dynamics are not known, a more general* 
and interesting case since it demonstrates the dependence on soil type, 
we have, from combining equations (3-4) and (3-5) , 


R 


£ 

r 


Combining equations 
the developed pull 
2 .. 


J F a = (c aea' M ez tan W ,i(id!; (3 ' 7) 

(3-2) and (3-7) results in a general expression for 



i 

(c a es w ez tan 4 es )d9,iz 


- w 


2 


(P+ma) 


2 


(3-8) 
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which is valid for any pressure distribution and any soil obeying 
the Mohr-Coulomb failure criterion. 

From this point the derivation will be continued for a cylindrical 
wheel, and it will be assumed that the pressure distribution is a representa- 
tive average with respect to z. Equation (3-8) then reduces to 

(P-taa) 2 - [^/(c ae .'*o e tan« to )d9 2 - W 2 (3-9) 

It now remains to evaluate the force components associated with R, 
and r£. 

FORCE COMPONENT EQUATIONS 

Consider the wheel-soil contact shown in Fig. 3-3 and the forces 
acting there. 

Depending on the sign of 9 (the location on the wheel-soil contact 

of the point in question) , the sign of the dN y forces will always be 

positive, while the sign of the dF , dF , and dN forces will vary 

y x x 

depending on 9 and the state of slip of the wheel. By appropriately 
adding up all the x and y force components, 

R 2 = (F + N ) 2 + (F + N r ) 2 (3-10) 

x x y y 

the vectorial resultant. 

The terms in equation (3-10) can be evaluated from: 

dN = OgdA = r0 Q d0 (3-11) 

* 

dF = TgdA - rTgd0 (3-12) 

} 
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F x = r j TgCos0d8 


N 


■ - r /°6 


sin6d0 


Fy=r TgSin0d0 


(3-13) 


N 


r 

-rj a ( 


y l 9 
SLIP AND INSTANTANEOUS WHEEL VELOCITY 


cos0d0 


The instantaneous velocity of the wheel, can be easily determined 

(see also Andreyer, Sitkei, and Janosi (Bekker, 1969)). 

Velocity components of a rolling wheel are shown in Fig. 3-4. 

The wheel is moving to the right with a velocity, v c / % » of the 

center of the wheel with respect to an instantaneous coordinate system, 

x, fixed to a far point in the soil. Point, p, has a tangential 

velocity with respect to the center of the wheel of v The desired 

velocity, V p/ X > * s that of the point p with respect to x and it is the 

vector sum of v , and v , . 

p/c c/x 


V p/x V p/c + V c/x 


(3-14) 


or. 


v p/x = V (v c/X - v p/c cos0)2 + (v p/c Sin6)2 (3 " 15) 


Slip may be defined as 


i = 


DR - DT 
DR 


(3-16) 


where DR = distance wheel would have revolved had it been rolling on a hard 
surface without slip, and DT = distance wheel actually travelled. Then 


DT 

DR 


« (1 - i) 


(3-17) 
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but also. 


PI V C/X 1 \U 

DR v . t v , 

p/c p/c 

where t = time; therefore, 


v c/x ■ v p/c (1 ' 1) 


Substituting equation (3-19) into equation (3-15) gives 


(3-18) 


(3-19) 


T p/x ■ ''p/cV» 1 - i > - cos9]2 + sln * 9 


(3-20) 


which defines the magnitude of the instantaneous velocity of any point 
on the wheel perifery. 

The direction of v , can also be determined. First, let 
p/x 




(3-21) 


, /v , - v . COS0 \ 

— ) 


V p/c 


(3-22) 



(3-23) 


Substituting equation (3-19) into equation (3-23), and simplifying gives 
. r> „ “I ( (1-i) - COS0 /o o/x 

-A X = 6 - tan [ ) (3-24) 

To establish a graphical solution for the instantaneous center of 
rotation, it would be interesting to determine at what vertical distance, 
y, from the bottom dead center of the wheel a line drawn at an angle A^ 
to the tangent at a point would intersect a vertical drawn through C. 
Referring again to Fig. 3-4, 
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but 


y " d 


d = 


r s in0 tan ( 0+A ^ ) 


= r sin0tanA o 


a = tan -l f (l-i) - cos 9 . 

A 2 tan ' 5H8 } 


z = r(l - cos9) 


With these substitutions 


(3-25) 


y = rsinStan [tan~ 1 (- ( -- ~ ■ -^"q- 036 ) ] - r(l-cos0) (3-26) 

Simplifying 

- sin0 (— — C0S Q) - 1 + cos0 

r ■ sin0 

or J - -i (3-27) 

This equation shows that the instantaneous center of rotation is located 
a vertical distance, ri, from the bottom dead center of the wheel-soil 
contact. For positive slip (+i) , a distance +y = +ri above, and for 
negative slip (-i) , a distance -y = -ri below the bottom dead center of 
the wheel. 

Therefore, the magnitude of the velocity of any point along the wheel 
perifery with respect to a fixed coordinate system in the soil can be deter- 
mined from equation (3-20) , and the direction of the same velocity can be 
determined using equation (3-24) or graphically using the instantaneous 
center of rotation. 


THE FRICTION CIRCLE METHOD 

The forces acting along a sliding contact are illustrated in 

t 


Fig. 3-5. 
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Clearly the direction of dR for any differential contact area, dA, is 
defined by the angle, ® 


tan® 


dF 

dN 


T 0 dA 

o e dA 


( 77 ^ + tanfi) 
a 0 


(3-28) 


Then, the moment arm of any dR (Fig. 3-1) must be 


r f q = rsin® (3-29) 

Noting that if the ratio C a /°Q in equation (3-28) is insignificant or 
zero, and if we have soil-to-soil sliding with friction angle <|> rather 
than 6 , 

r^g - rsin® ■ rsin(tan ^tan<|>) » rsinij) = r^ (3-30) 


Equation (3-30) defines the radius of the friction circle as used by 
Taylor (1937) . (Taylor first introduced the friction circle method for the 
analysis of the stability of cohesionless embankments.) In equation (3-30), 

r^g is independent of 0 and therefore equal to r^. 

Combining equations (3-4) and (3-29) we obtain 


r^g sin® R 


(3-31) 


Equation (3-31) expresses the extent by which the moment arm to the 
resultant R differs from the moment arm to a dR force. That there is a 
difference between r£ and r^g can be clearly illustrated for c & = 0 as 
shown in Fig. 3-6. Whereas all the dR forces in Fig. 3-6 are tangent to 
the friction circle with radius r^g, the resultant of any two dR forces, 
such as dR^s acts at a slightly greater distance, rj^* Analogously, the 
resultant of all dR forces will also act at a somewhat greater distance, 
r^. The difference is anticipated to be small, however, and the correction 
factor, r^/r^g, will be evaluated for certain cases of interest. 
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RIGID WHEEL IN COHESIONLESS SOIL 


For a rigid cylindrical wheel driven at constant velocity in a 


cohesionless soil for which c =0, equation (3-9) reduces to 

SL 


p 2 


-B /• 


tan6 Qs ff Q d0 


- 


(3-32) 


and equation (3-31) becomes 



1 

sinS Q 

0s 


F 

a 

R 


(3-33) 


Using equations (3-5) and (3-33) , equation (3-32) can be expressed in 
terms of the r^/r^ ratio 


P 2 


rVr 

r f /r f 



- WT 


(3-34) 


Using the force component equations (3-10) , (3-12) and (3-13) , the 
correction ratio, r^/r f , can be solved from equation (3-33) for any 
integrable distribution of radial pressure at the wheel-soil contact. 
For a constant pressure distribution, for example, equation (3-33) 
reduces to 



0 


T 


'll /l-cosQ^ 


(3-35) 


which is identical to Taylor's solution (0^ = total wheel-soil contact 
angle) . Taylor (1937) presents values of this correction ratio (Fig. 3-7) 
for both a constant and a sinusoidal pressure distribution. A constant 
pressure distribution is a good approximation for a pure clay; a sinusoidal 
pressure distribution is a good approximation for a pure sand. Other 

i 

pressure distributions would likely lead to some intermediate correction. 
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Since Fig. 3-7 shows that the correction factor is very small (1-6%) 
for either case, for contact angles less than 70°, we may consider 
equation (3-34) a nearly rigorous solution for the pull that can be 
developed in a cohesionless soil. 

The theoretical developments presented in this chapter are all 
based on the Mohr-Coulomb failure criterion, and depend on a knowledge 
of the shear and normal stress distributions along the wheel-soil 
contact. It is recognized that these pressure distributions are 
difficult to determine. In the next chapter it will be shown that, for 
all practical purposes, certain observations make it possible to bypass 
the difficult question of pressure distributions. 
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Fig. 3-7. Correction to be used with the 
friction circle method. 


(Taylor, 1937) 
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CHAPTER 4. APPROXIMATE WHEEL-SOIL INTERACTION THEORY 

LINE OF ACTION OF RESULTANT OF RADIAL STRESSES— A FUNDAMENTAL OBSERVATION 

The line of action of the resultant of radial stresses approximately 
bisects the wheel-soil contact angle, 0^, for all values of slip, s. 

This observation is supported by: 

a) A rational argument based on a cylindrical wheel rolling down 
a soil slope at constant velocity 

b) An empirical argument based on data from spheres rolling down a 
soil slope at constant velocity 

c) Test data from many wheel-soil interaction experiments. 

A relatively simple approximate theory is then developed. 

EVIDENCE IN SUPPORT OF THIS OBSERVATION 
Cylindrical Wheel Rolling Down a Soil Slope 

A cylindrical wheel rolling down a soil slope is shown in Fig. 4-1. 
The symbols in Fig. 4-1 are defined as follows: 

W = weight of wheel 
m = mass of wheel 
R = resultant soil reaction force 
R = component of R parallel to slope 

X 

Ry = component of R normal to slope 
r = radius of wheel 

r^ = distance from center of wheel to line of action of R 
I = moment of inertia of wheel 
a = linear acceleration of wheel 
= angular acceleration of wheel ' 


u 
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x,y = coordinate directions 
a = slope angle 

Negative slip usually takes place as a wheel rolls freely down a slope, 
and, therefore, a / ur. 

Applying Newton's Second Law, using the method of dynamic equilibrium, 
ma and Iu are considered acting opposite to their actual sense, as shown by 
the dotted arrows in Fig. 4-1. We then have 


++ 


R = W cos a 

y 


(4-1) 


E R = 0 = R - W cos a, 

y y 

+ V R = 0 = W sin a - R - ma, R = sin a - W^- (4-2) 

Z_* x x x g 


-hi Y, M g = 0 = Rr f " Iu > 


The resultant, R, in terms of R and R , is 

y x 

R 2 ■ R 2 + R 2 


R =|W^ 
5 r f g 


(4-3) 


(4-4) 


Substituting R and R from equations (4-1) and (4-2) into equation (4-4) 

y x 


gives 


R = W £^3in a - + cos 2 aj 


1/2 


(4-5) 


If a and u were established experimentally (this can be easily done by 
taking movies of the rolling wheels), R and r^ could be determined from 
equations (4-3) and (4-5) . This would give the magnitude and direction 
of the resultant soil reaction, R, as well as its point of action on the 
wheel surface. 

For constant velocity rolling, (a = 0), equation (4-5) shows that 
R = W, and has the same line of action. Net shear stresses at the wheel- 
soil contact must be zero since moments about the wheel center are zero. 
Because soil contact pressure tends to distribute itself about the point 
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of application of the load, W, the resultant of soil stresses, R, must 
pass approximately through the middle of the wheel-soil contact area. 

Further, the continuous process of rolling requires that the line 
between the wheel center and the wheel-soil rear contact be nearly perpendi- 
cular to the slope. (Soil exhibits some rebound, and the line between the 
wheel center and the wheel-soil rear contact deviates from a perpendicular 
by B^). There conditions are illustrated in Fig. 4-2. From these condi- 
tions a unique relationship can be established between track geometry and 
slope angle. If the line between the front and the rear wheel-soil contacts 
is horizontal, 

f = 1 - cos (2a + e 2 ) (4-6) 

It has been argued that for a cylindrical wheel rolling at constant 
velocity down a slope, the resultant of radial contact stresses bisects 
the wheel-soil contact angle, 0^, and a unique relationship exists between 
sinkage and slope angle (equation (4-6)), which is independent of wheel 
load. A free-body diagram of a wheel rolling on a horizontal surface 
would differ from the above case only with respect to the direction of W, 
resulting in R = W/cosa. 

Therefore, a rational argument exists for the validity of the observa- 
tion stated on page 29 . What about experimental evidence? 

Spheres Rolling Down a Soil Slope 

For a sphere (spherical wheel) rolling down a soil slope, the track 
width, w, is related to the sinkage geometry (Fig. 4r2) by 

sin(2a + 0„) = = — 

2 r D 


(4-7) 
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which, again, is a unique relationship independent of sphere load, for 
constant velocity rolling only. 

Spheres rolled on Yuma sand provided an opportunity to test equation 
(4-7) for ©2 = 0. These experiments are described in detail by Hovland 
(1970). Fig. 4-3 shows the comparison. The agreement is best for loose, 
moist sand. In this material the tracks left by the rolling sphere had 
sharp, well-defined edges which could be accurately measured. The devia- 
tion is greatest for loose, dry sand where the tracks tended to cave or 
slump, and where the crest-to-crest track width, w, measured was undoubtedly 
somewhat greater than the track width at the time the sphere was in contact 
with the soil. 


Experimental Data on Wheel-Soil Contact Angles 
and Radial Stresses 

Physical reasoning would indicate that the resultant of radial stresses, 
at least for low values of slip, should approximately bisect the contact 
angle. The value of stress must be zero at the front contact and zero at 
the rear contact. For high positive or negative values of slip, however, 
this is not obvious. To compare data from various experiments for all 
values of slip, a redefinition of slip is first desirable. 

The most common definition of slip, equation (3-16), is 


i 


DR - DT 
DR 



(3-16) 


where DR ■ distance wheel perifery rolled, and DT ■ distance wheel traveled. 
Fig. 4-4a shows a plot of equation (3-16). A purely skidding wheel, DT/DR 
= 00 , would plot at i = - 00 . This is consistent with the instantaneous 
center of rotation, but not useful for developing theory valid for the 

whole slip range. Some investigators (for example, Sela, 1964) define 

/ 
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0 10 20 30 40 50 

Slope Angle, a, degrees 


Fig. 4-3. w/D vs. slope angle, a, for 
constant velocity rolling. 
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positive slip by equation (3-16) but negative slip by 


• - PT ~ DR _ DR 
-1 “ DT “ 1 ~ DT 


(4-8) 


Therefore', it is necessary to have one mathematical expression which 
is valid for all values of slip. This is achieved by 

i DT 

DR - DT ” DR ,, „„ 

s = = — (4-9) 


DR + DT 


1 + 


DT 

DR 


The author prefers equation (4-9) because a purely spinning wheel has slip 
s = +1, and a purely skidding wheel has slip s = -1 (see also Fig. 4-4b) . 
Since all definitions are functions of DT/DR, one can always go from one 
definition to another. Occasionally in the literature the definition for 
slip has not been explicitly stated. In such a case, it has been assumed 
that slip was defined by equation (3-16) . 

Of the data considered, Sela's (1964) showed the most consistent 
relationship and also extended over a larger range of slip values (Fig. 
4-5) . The data suggests that the position of the resultant of radial 
(normal) pressures, 0^, moves slightly forward with increasing slip. All 
the data (Fig. 4-6) support the same trend. 

A statistical analysis of all the points between s = - .3 and s = +.3 
(Fig. 4-7) shows that the mean is 0.49 and the standard deviation is 0.06. 
Comparing information in Figs. 4-5, 4-6, and 4-7, an average drawn through 
Sela's data is about one standard deviation above the mean based on all 
the data, but the slope of the average line is about the same. 

In order to use the data in Fig. 4-6, for example, it is necessary 
to separate out 0^ and 02* Considering © 2 » when s = -1, 02 should equal 
zero, and when s = +1, ©2 should approach 0^ as slip sinkage increases. 
Thus, if were Pl° tte d vs. slip, some kind of curve beginning with a 



Sela, 1964 

rigid wheels and dry sand 
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Fig. 4-5. Line of action of the resultant of radial stresses as a function 
of slip (Sela's, 1964 data). 


Uffelmann, 1961 clay rig. wheel o Sela, 1964 sand rig. wheel 

Sohne and Sonnen, 1961 silty sand rig. wheel O Onafeko and Reece, 1967 sand rig. wheel 

Vincent, 1961 sand pneu. tire • Krick, 1969 sandy loam rig. and pneu. 



f 


Fig. 4-6. Line of action of the resultant of radial stresses as a function of slip 
(all collected data) . 



Number of Data Points, -.3<s < +.3 
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low value of t ^ e ne g at ive slip range, and increasing with increasing 

slip should be expected. Fig. 4-8 shows an attempt to plot vs ‘ 

slip. The scatter does not invalidate the plot, but it does illustrate 
that much more detailed experiments are necessary to establish what 
influences 0^5 particular attention should be given to the influence of 
soil types. 

Fig. 4-9 shows a plot of 0 /0. vs. slip. The plot in Fig. 4-6 shows 

N i 

a better correlation, however. 

The data in Fig. 4-6 can be approximated by the line 


0 N + 9 2 


9 N + 9 2 
9 1 + 9 2 


- 0.5 + 0.1s 


(4-10) 


A curve through the centroid of points in Fig. 4-8 can be approximated by 


s = -0.5 + 40 


e 2 .e 


(4-11) 


or 


^2 _ e 2 . 3 + .5 .1/e 

0 T ' 6 1 + 0 2 40 

Combining equations (4-10) and (4-12) 


(4-12) 




- (.5 + .Is) + 


(.5 - .Is) 


1 - ( 


40 

s+.5' 


1/e 


(4-13) 


Judging by the fit between the data and the curve (equation (4-13)) in 
Fig. 4-9, the equations appear to be adequate, and the curve in Fig. 
4-8 is a good first approximation. 


THEORY 

Forces Acting on a Wheel 

The forces acting on a wheel operating on a horizontal surface are 



Uffelmann, 1961 clay rig. wheel o Sela, 1964 sand rig. wheel 

Sohne and Sonnen, 1961 silty sand rig. wheel O Onafeko and Reece, 1967 sand rig. wheel 

Vincent, 1961 sand pneu. tire ® Krick, 1969 sandy loam rig. andpneu. 



Fig. 4-8. Rear wheel-soil contact angle as a function of slip. 




Uffelmann, 1961 ciay rig, wheel o Sela, 1964 sand rig. wheel 

Sohne and Sonnen, 1961 silty sand rig. wheel O Onafeko and Reece, 1967 sand rig. wheel 

Vincent, 1961 sand pneu. tire • Krick, 1969 sandy loam rig and pneu. 


43 



Fig. 4-9. The ratio of to front wheel-soil contact angle vs. slip. 
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shown in Fig. 4-10. Four cases are shown: 1) a driven wheel generating 

pull (+F, +P) , 2) a self propelled wheel (+F, P = 0) , 3) a towed wheel 
(F = 0, - P) , and 4) a braked wheel (-F, = - P) . As shown in Fig. 4-10, 
if 0^ and W are known and F is assumed to act perpendicular to N, the 
problem is determinate. Arguments and data were presented for 0^ in 
the previous section. F does not always act perpendicular to N, and a 
correction will be developed. Interestingly, the line of action of F is 
likely to deviate most from a perpendicular to N for low and moderate 
mobilizations of F. For high values of ±F, when consequences of the devia- 
tion could be significant, F acts approximately perpendicular to N. 

The following expressions can be written directly from Fig. 4-10: 

(4-14) 
(4-15) 
(4-16) 

In these equations, all forces are vectorial unless otherwise noted. 

From equations (4-14) and (4-15) , 

P 2 = N 2 + F 2 - W 2 (4-17) 

The algebraic sum of shear stresses at the wheel-soil contact can be 
expressed by 

F a = Ntan6 0s + Ac a0s * (4 - 18) 

where A = the wheel-soil contact area = wQ^r (w = wheel width) . (To 

abreviate slightly the subsequent equations, the subscripts 0 and s are 

omitted from 6 and c ) . 

a 

Taking moments about the center of the wheel, the v^corial F is 

F 

a 


R 2 = N 2 + F 2 


R 2 = P 2 + W 2 


N - Wcos0„ - Psin9„ 
N N 


r/r 


The correction ratio r/r will be considered separately. 


(4-19) 
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Substituting equations (4-16), (4-18), and (4-19) into equation (4-17) 

gives 



For a pure sand (c = 0) , equation (4-20) becomes 

Si 



For a pure clay (6 = 0) , equation (4-20) becomes 

(I)' ' ( c ° s6 N - I + {—) - 1 


(4-21) 


(4-22) 


For a towed wheel, net shear stresses are zero, which caft also be accom- 
plished by letting c & and 5 in equation (4-20) equal zero 



(4-23) 


These equations are valid for all conditions illustrated in Fig. 4-10. 
Although the equations are lengthy, they are not difficult to solve; the 
variables are 0^, 6, c^, and 0,^. Testing yields 6, and c a , for every case 
of interest. Anticipated sinkage and previous developments in this report 
yield 0^ and 0 T » 


Consideration of Magnitude and Direction of F 

It was noted by equation (4-18) that F q in that equation was algebraic, 
and that, taking moments about the center of the wheel, the vectorial F is 


F = 


r/r 


> 


(4-19) 
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where r 


moment arm to action of F. Then 
F 


r_ 

r 


/ F 2 + F 2 
x y 


(4-24) 


or 


r T 0 d0 


r y (r / T 0 cos0d0) 2 + (r / T fl sin0d0) 


(4-25) 


0" 


For a constant shear stress distribution, equation (4-25) reduces to 

0 _ 


r 

r 


/l - 2 cos 


(4-26) 


which is identical to equation (3-37) . Therefore, the same correction 
chart (Fig. 3-7) can be used. 

For a sinusoidal shear stress distribution, equation (4-25) reduces 


to 


- !-£)' 


e„ 


(4-27) 


cos 


which is, again, the same correction obtained by Taylor (1937) the 
correction is found from Fig. 3-7. 

The amount by which the direction of F deviates from a perpendicular 
to the direction of N is illustrated in Fig. 4-11. This deviation can be 
expressed by 

F 1 

(4-28) 


F 1 

l 0 N + 9 2 " 9 f! = l 0 N + 0 2 " tan " 

x 


where 


F * r/t Q sin0d0 
_JL_ = 9 

F 1 r /x o cos0d0 
x o 


(4-29) 


Solving equation (4-28) , using equation (4-29), for a constant and a sinu- 
soidal shear stress distribution respectively, gives 
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Tq = constant, 


Tg = sinusoidal. 



(4-30) 


In other words, for constant and sinusoidal shear stress distributions, 
equations (4-30) show that the deviation is equal to the deviation of the 
line of action of N from 0^/2. Figs. 4-6 and 4-7 show that 0^ + 0^ 
is not likely to deviate much from 0 ^/ 2 . 

Inaccuracy Resulting from an Error in Selected Contact Angles 

Observing equations (4-20) through (4-23) , it appears that an error in 
selecting 0 ^ and associated ©2 would have the greatest influence in the 
case of a towed wheel, equation (4-23). Therefore, only the case of a 
towed wheel will be considered. This case can be expressed more simply 
from Fig* 4-10 as 

**T 

-y- » tan0 N (4-3D 


The magnitude by which the computed pull differs from the correct pull 
can be expressed by the ratio 


(P T /W) 

(P T /W) 


comp. 

correct 


tan( 0 M ) 

N comp. 

tan( 0 M ~) 

N correct 


(4-32) 


Based on the data in Figs. 4-6 and 4-8, it is most likely that a high 
value of (0 N + ®2^®T * S assoc * ate ^ with a high value of 62 / 0 ^. Therefore, 
the deviation of 0 ^ from the mean can be expressed as a function of 0 ^ 


®N + 0 2 


'N 


- mean ± standard deviation 


and 02 from 
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(4-33) 


.56 - .3 = .26 
.50 - .2 = .30 
.44 - .1 = .34 

The pull ratio computed from equation (4-32) , using the mean 
(0 = 0.3 0,j,) as the correct value, is presented in Table 4-1. 

Table 4-1 Ratio of Computed to Correct Pull 


9 N +e 2 

0 T 



0 

T 



0 

20° 

O 

O 

60° 

90° 

120° 

high 

0.5+0.06 

1.00 

0.87 

0.86 

0.86 

0.85 

0.84 

low 

0.5-0.06 

1.00 

1.13 

1.14 

1.15 

1.16 

1.18 


Since for a normal distribution, the range from -(one standard devia- 
tion) to +(one standard deviation) contains 68.27% of the data, it can be 

f 

concluded that for about 68% of wheel-soil interaction cases the computed 


O'- ' —‘O 


0 N 0 2 

A— = (.50 ± .06) - —■ 

0 T 0 T 


Then we have for 


high 


0 N + 0 2 9 N 

e„ * e m 


mean 


N 


low 


N 


’ 9„ 
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pull will be within about 15% of the correct value. This conclusion is, 
of course, based on the data presented in Figs. 4-5 through 4-9. With 
future careful testing of (0^ + ^^^T aS a ^ unct ^ on s °il conditions 

and wheel conditions, it should be possible to reduce this theoretical 
error to a smaller value. 
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CHAPTER 5. WHEEL SINKAGE 
WHEEL SINKAGE AND CONTACT ANGLES 

The towed force (negative pull = -P = P T > expressed by equation (4-23) 

is a function of W and 0^ only. As illustrated in Fig. 4-10, every 

state of slip of a wheel has associated with it a motion resistance, N x< 

For a horizontal surface the N /N ratio can be expressed directly 

x y 


= tan0, T (4-31) 

The sinkage, z, associated with this angle, 0^ , is illustrated in 
Fig. 5-1. Sinkage is here defined as the vertical distance between the 
bottom dead center of the wheel and the front wheel-soil contact 

z = r(l - cos0^) (5-1) 

from which 

0^^ = cos" l (l - 2 |) (5-2) 

where D = wheel diameter. From Fig. 4-9, we note that 0^ /0^ is 
nearly constant and equal to 0.4 for most likely slip values. Therefore, 

0 N =0.4 cos"" 1 (1 - 2 |) (5-3) 

Combining equations (4-31) and (5-3) gives 

^ = tan [0.4 cos-^l - 2 ^) ] (5-4) 

y 

Equation (5-4) relates sinkage to the motion resistance for a rigid 

cylindrical wheel. It is an approximate solution essentially independent 

} 

of slip. Improvements to equation (5-4) should be sought from improved 
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data of the type plotted in Fig. 4-9. To check equation (5-4), experimental 
data from towed wheels are compared in Fig. 5-2, 

Many relationships between sinkage and rolling resistance have been 
advanced previously, and many of these are conveniently summarized by 
Schuring (1972). He defines rolling resistance as p = sin0^. With this 
definition and equation (5-3), rolling resistance can be expressed as a 
function of sinkage. This leads to results which are almost identical with 
the Bernstein line (Schuring, 1972). The experimental data reported by 
Schuring in this paper deviate from the Bernstein line particularly for 
driven, pneumatic tires and low sinkage. (Actual rolling resistance appears 
to be considerably less than predicted.) Schuring proposes an explanation 
for this deviation: For low sinkage, the rear angle, tends to be 

proportionately larger, which leads to a smaller 0^ and smaller rolling 
resistance. He also proposes a relationship which seems to account for the 
influence of 02 rather well. 

PRESSURE-SINKAGE RELATIONS 

While equation (5-4) relates sinkage to motion resistance, N^, it 
cannot be used to compute sinkage directly from soil parameters. To do 
this, an independent relationship is required. 

Perhaps the most commonly used formula In the field of off-road 
mobility, to define vertical stress-strain relationships, is (Bekker, 

1969) 



where q = bearing pressure on a plate of width, w, and sinkage, z. The 
constants k c and k^, and exponent n are determined from testing and curve 
fitting. Bekker (1969) presents a thorough description of the use of 
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z/D 


Fig. 5-2 Ratio of towed force to wheel load vs. ratio of sinkage 
to wheel diameter. 
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equation (5-5), and many experimentally determined pressure-sinkage curves. 

A significant amplification of the type of pressure-sinkage relation- 
ship represented by equation (5-5) has been made by Firth (1968) . On the 
basis of dimensional analysis, he derived a more general pressure-sinkage 
relationship . 



where V = volume of probe below plane of original surface 

A = cross-sectional area of probe in plane of original surface 
Jl = major dimension of area, A 

m and n are exponents to be determined from testing 
Some separately postulated relationships, such as equation (5-5), are, in 
fact, special cases of equation (5-6) (Firth, 1967). 

Alternatively, it is possible to use plasticity theory for pressure- 
sinkage relationships. A notable application of plasticity theory to the 
wheel-soil interaction problem has been made by Karafiath (1971) in an 
effort to predict stress distribution beneath wheels. His studies have 
also provided valuable information on the effect of soil pore water 
pressures in wheel-soil interaction (Karafiath, 1972). 

Bearing capacity theory (a special solution based. on plasticity theory) 
can also be used for pressure-sinkage relationships. The bearing capacity 
equation 

Y w 

q = ~T~ N y + cN c + q ' Nq (5-7) 

* 

was first derived by Terzaghi (1943) . Equation (5-7) is a function of 

sinkage by 1) q’ = y z/2, and 2) q = N/A, where A, the wheel-soil contact 

s 

area, is a function of sinkage. 

f 
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The main assumption in both Karafiath’s analysis and equation C5--7) 
is that of two-dimensionality (an infinite dimension perpendicular to 
the plane of rolling in the case of wheel-soil interaction). Analysis 
based on this assumption, therefore, may be, at best, adequate for 
relatively wide wheels and for those cases where theory is substantiated 
or modified by experiments. 

It is customary to modify the bearing capacity equation by shape 
factors (s y and = s^) for other than an infinitely long load-soil 
contact area. Using shape factors, and noting that = N c tan<])+-l, 
equation (5-7) can be expressed for rectangular loaded areas as 

y w 

q = -§- s y N y + s cq [(c+q'tan<!>)N c + q') (5-8) 

where y g = soil unit weight, w = width of contact area, q’ = surcharge, 
and N y , N c , and are bearing capacity factors which are functions of 
soil friction angle, cf>, and slope angle, a, (Meyerhof, 1951). 

While equation (5-8) has some empirical backing and considerable 
precedent in soil mechanics practice, it has not, to the author’s 
knowledge, been experimentally modified or substantiated for rectangular 
wheel-soil contact areas. Although equation (5-8) will be suggested 
for long, rectangular, wheel-soil contact areas, its use in wheel-soil 
interaction must be considered hypothetical. 

Bearing capacity theory can, however, be modified by experiments 

for specific cases of wheel-soil interaction. Extensive testing (Hovland, 

1970) demonstrates that bearing capacity theory can be used to evaluate 

pressure-sinkage relationships of rolling spheres and spherical wheels. 

Fig. 5-3 shows a comparison between theoretical (modified bearing 

capacity theory) and experimental values of the density ratio , y /y , 

r s 



a 25° 
o 35° 


4>= 37° 
X s = 1.465 
c = 0 


+ o 



0.2 0.4 0.6 0.8 0 
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(Loose air -dry sand) 


$ 3 42° 
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w/D 
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(Dense air- dry sand) 


•+ 



<t> = 375® 

7s 3 1.47 

C * 1.07 KN/m 2 
w * 1.50% 
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w/D 

(c) Car 4 

(Dense moist sand) 


/ <f> *32° 

/ + y s = 1.25 
/ c = 0.46sct^rf 
w =1.55% 

0.2 0.4 0.6 0.8 1.0 

w/D 
(d) Car 5 

(Loose moist sand) 


Density ratio vs. w/D ratio for spheres rolling 
on Yuma sand (curves are theoretical for constant 
velocity rolling; points are experimental) . 
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vs. track width to diameter ratio, w/D. Fig. 5-4 shows a comparison 
between theoretical and experimental values of sphere diameter vs, track 
width to diameter ratio, w/D. In these figures = sphere unit weight. 
Fig. 5-5 shows a comparison between theoretical (calculated) and ex- 
perimental friction angle values. The theoretical values were again 
calculated from bearing capacity theory modified to the rolling sphere- 
soil interaction problem (see equation (7-12)). 

Such testing showed that, for the sphere-soil interaction problem, 

the bearing capacity factors (Meyerhof, 1951) need to be modified in 

equation (5-8). Specifically, for a sphere or spherical wheel 

.. -IN 

N = — v 

ysphere 4 

and N = - N 

c , 2 c 

sphere 

Substituting these bearing capacity factors into equation (5-8) gives 
for a spherical wheel 

Y w N r N -i 

9 = -f“ s y X + S cq L (C + q ' tan & T + q J (5_9) 

We have, therefore, matched bearing capacity theory in an approxi- 
mate way to two end-points of the wheel-soil interaction problem; 1) 
a rather narrow cylindrical wheel for which equation (5-8) applies,^ and 
2) a spherical wheel, which is analogous to a wide wheel or a very light 
wheel, for which equation (5-9) applies. What about intermediate 
conditions? 

Fig. 5-6a illustrates wheel-soil contact areas, as viewed from above, 
for various wheels : 

i 


(1) A spherical wheel. 






Calculated Friction Angle 


o 

o 


No. of data points 


/\ 45 ° 

1 1 

20 30 40 50 

Triaxial Test Friction Angle 

(a) All Data 



Triaxial Test Friction Angle 
(b) Constant Velocity Data 


Fig. 5-5. Comparison between calculated and triaxial test 
friction angle values. 




Direction of motion 


(I) (2) (3) (4) (5) 



a) Wheel-soil contact areas for various wheels. 



b) Correction factor to the bearing capacity factors of Eqs. (5-8) 
and (5-9) . 


Fig. 5-6. Suggested Correction of bearing capacity theory for 
wheel-soil interaction. 
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(2) A wide cylindrical wheel with a contact area approximately 
equal in shape to that of a spherical wheel. 

(3) A cylindrical wheel with a square contact area, at which 
point the Jl/w ratio in the shape factor expressions 


and 


Y 


1 - 




S = 1 + 

cq 


3 


~) 

w i 



(5-10) 


(4) 


equals one. 

A typical cylindrical wheel with 
£/w ratio to be larger than one, 
must be inverted, giving 

and 


sufficient sinkage for the 
at which point the ratio 


(5-11) 


(5) A rather narrow cylindrical wheel. 

For the condition of Fig. 5-6a (1), equation (5-9) applies, using 
the ratio £/w in determining the shape factors. For the condition of 
Fig. 5-6a (5), equation (5-8) applies, using the ratio w/Jl in determining 
the shape factors. For intermediate cases, compute the shape factors 
from equations (5-10) and (5-11), and select from Fig. 5-6b the correction 
factor to the bearing capacity factors of equations (5-8) and (5-9), In 
all cases, select the bearing capacity factors from Meyerhof (1951) 

(Figs. 5-7, and 5-8), using 0^ = Meyerhof's 8 m * 
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CHAPTER 6. MOBILIZATION OF SHEAR STRESSES 


SLIP AT A POINT 

The relationship between measured slip and slip at any point along 
the wheel-soil contact is illustrated in Fig. 6-1. Slip is established 
from measured distance through which the wheel revolved, DR, and distance 
through which the wheel traveled, DT. The measured slip (slip by its usual 
definition) takes place at the bottom dead center of the wheel, and may 
be expressed as 


DR - DT 
DR + DT 


(4-9) 


from which 


DT _ 1 - s 
DR " 1 + s 

Slip at a point, s^, is a function of 0, and may be expressed from 
Fig. 6-1 as 

s - DR - DT cos8 
0 = DR + DT cos0 

combining equations (6-1) and (6-2), 

- 1 ~ COS0 + s (1+CQS0) 


( 6 - 1 ) 


( 6 - 2 ) 


(6-3) 


0 1 + COS0 + s (1-COS0) 

Because the expression of slip at a point (equation 6-3) may be new, 
it is desirable to see if predictions of its value are realistic. A 
surface of slip at a point, Sq, is plotted in Fig. 6-2. The surface 
appears to be realistic and in accord with physical reasoning. Note 
particularly the cross-hatched area at the towed point; as the wheel load 
increases, the wheel-soil contact angle increases, and the towed point 
moves to a higher negative slip value. When the necessary relationships 
are established, later in this chapter, the towed point will be reconsidered. 
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Fig. 6-2. Surface of slip, Sq, at any point. 
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SOIL DISPLACEMENT VS. SLIP 

A shear ring type of test (Bekker, 1969) (Fig. 6-3b) is probably most 
realistic for establishing wheel-soil interface strength parameters c 

3 . 

and 6. In this type of test, soil is allowed to strain 

analogously to what happens to the soil under a wheel. Compare this test for 
example, to a direct shear test where strain is limited to a very narrow zone; 
peak strength would occur at a much smaller displacement in a direct shear 
test than in a ring shear type of test shown in Fig. 6-3b. 

In this type of test, however, as in the case of the wheel, it is 
difficult to establish fundamental stress-strain relationships since 
strain cannot be defined. That is, the vertical distance, Az, to which 
soil deforms is unknown. It is, however, possible to establish displace- 
ment relationships which may be useful. 

Fig. 6-3 illustrates soil displacements and terminology for 
wheel- and ring-soil interaction. For the wheel, the rate of displacement 
can be expressed as 


A 0 r 9 r (1-i) 0 r 

w DR-DT ww-ww ww 


(6-4) 


w 


w 


w 


w 


where t = time increment for wheel, and i = slip as defined by equation 
(3-16). 

For the ring shear test (Fig. 6-3b) 


\ 9 t r t 
c t ' c t 

where t = time increment for test. 

t * 

An equivalence may be written from equations (6-4) and (6-5) by 
assuming that the rate of soil displacement for the wheel and the test 
are the same 


(6-5) 
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b) Shear Ring Test (three dimensional view) 


Fig. 6-3. Soil displacements for wheel and test. 
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A /t 0 r t 

w ” = 1 - “ . r (1) 




0 * t 

t t w 


( 6 - 6 ) 


If Instead it is required that the strains for the wheel and the test be 
equal. 


e 9 , r Az 

JL m i = JL_ w t .. 

e . 9„ r„ Az w 

t t t w 


(6-7) 


While these equations may be useful in selecting test setup and 
equipment, they do not tell us the slip, associated with a particular 
displacement. 

From equation (6-4) 

A - 6 ri (6-8) 

w w 

where r is the radius of the wheel, consistent with our earlier 
definition. Using equations (3-16) and (4-9), equation (6-8) can be 
expressed in terms of "s' 1 


A = 


29 r 
w 


w 1 


(6-9) 


+ 1 


Solving for s from equation (6-3) and substituting for s into equation 
(6-9), 


A = 0 r 
w w 


1 - 


g - s 9> 
(1+Sq) cos0 


( 6 - 10 ) 


Note that if 9 = 0, Sg = s, and equation (6-10) reduces to equation (6-9). 

To facilitate certain comparisons and illustrative calculations to 
be presented in this report, equation (6-10) is plotted in Fig. 6-4 for 
0=0, assuming 0^ = 0^. In reality, could have any value; however, 
physical reasoning suggests that 0 w = 0^ is a good first approximation, and 
it provides a convenient bridge from slip atp a point to displacement. 
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With better understanding of the wheel-soil interaction process, 
the functional relationship expressed by equation (6-10) undoubtedly 
will need to be modified. For instance, A is probably a function of 
soil type; a coefficient, which would be very small for a brittle soil 
and relatively large for a very plastic soil, may be required. Further, 
it appears that soil elements influenced by the wheel are subjected to 
various states of shear; while some soil elements are in a state 
analogous to direct shear for high values of slip, some soil elements 
are in a state analogous to compression for low values of slip. These 
complexities also need to be considered. 

WHEEL-SOIL INTERFACE STRENGTH PARAMETERS 

As in many other engineering problems, the fundamental relationships 
to be obtained from testing are stress-displacement curves. Suppose that 
testing of a wheel surface material to soil, using a ring shear device, 
yields the curves shown in Fig. 6-5a (three tests - three different normal 
pressures - are adequate). From Fig. 6-5a shear stress vs. normal stress 
lines can be plotted for various displacements (Fig. 6-5b) . From Fig. 6-5b 
strength parameters, c & and 6, can be plotted as a function of displacement 
(Fig. 6-5c) . These are the required wheel-soil strength parameters. 

SHEAR STRESS SURFACE 

Experimental wheel-soil interaction data suggest a shear stress 
surface, T = f (s, 8). An attempt will now be made to predict such a 
surface from the previously presented relationships. The procedure is 
as follows: 

1. Testing yields soil strength parameters, c and 6 as a function 

Si 

of displacement (Fig. 6-5) . / 
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2. Slip, s, can be converted to slip, Sg, using equation (6-3), 

3. Displacement can be related to slip, Sg, from Fig, 6-4. 

4. The shear parameters are then known for slip at any point, Sg, 
along the wheel-soil contact. 

5. Therefore, the shear stress mobilized at any point along the 
wheel-soil contact can be computed from 

T se - c as0 + °e tan6 se <6 " 11) 

In this prediction, the normal stress distribution, 0gi may be assumed 
sinusoidal for a sand and constant for a pure clay. For any real soil, 
the normal stress distribution is probably closer to sinusoidal. 

The predicted shear stress surface is shown in Fig. 6-6, which shows 
a three-dimensional view. In this figure the contact angle, 0, is plotted 
on the horizontal axis, slip, s, is plotted on the diagonal axis, and shear 
stress is plotted on the vertical axis. The front part of the surface 
(dotted lines) is below the 0, s plane and the rear part (solid lines) is 
above the 0, s plane. The normal stress distribution was assumed to be 
sinusoidal. Soil strength parameters were obtained from Fig. 6-5, with 
c = 10.0 at displacement * 0. In Fig. 6-4, 0 r was assumed to equal 

9 - 1 

10 inches. For a 20-inch diameter wheel, 0 T would then be approximately 
60 degrees. The rear wheel-soil contact angle, © 2 * was assumed to be 10 
degrees . 

The predicted shear stress surface appears quite realistic. The 
towed point for this case occurs approximately at s * -0.08; it is indicated 
by the cross-hatched area in Fig. 6-6. If friction at the wheel hub is 
zero, torque must be zero for a towed wheel (see also Yong and Webb, 1969), 
and the towed point is characterized by 


/ V 0 = 0 


(6-12) 
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In summary, given three stress-displacement curves, it is possible 
to estimate the shear stress mobilized at any point along the wheel-soil 
contact, for any value of slip. 
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CHAPTER 7. SOIL INERTIA IN WHEEL-SOIL INTERACTION 
In previous investigations of wheel-soil interaction phenomena, the 
influence of inertia forces of moving soil have usually been assumed to be 
negligible. Although this assumption is probably valid for many wheel- 
soil interaction problems, it is supported by limited investigations and 
experimental data. 

The investigations described in this chapter present a universally 
applicable hypothesis, to be tested experimentally. The hypothesis is 
based on the concept of a soil wedge being continuously formed and 
accelerated in front of the wheel. 

Fig. 2-1 illustrates a track and associated shear surfaces resulting 
from pulling a 24- inch (61-cm) diameter wheel on Yuma sand. The extent of 
the shear surfaces is evaluated, and theory is presented from which the 
mass and acceleration of the moving soil can be deduced, thus forming the 
basis for evaluation of soil inertia forces. 

THEORETICAL ANALYSIS 

The resistance to a rolling wheel caused by the inertia of the soil 

set into motion by the wheel, R^., is the product of the mass of soil 

involved, m , and its acceleration, a , 
s s 

R = m a (7-1) 

1 s s 

Mass of Moving Soil 
Boundary of soil wedge ; 

A mechanism of shear in front of and below a wheel rolling down a soil 
slope is shown in Fig. 7-1. Imagine that the wheel center is at position 1. 
As the wheel center moves from position 1 to 2 and then to 3, the front of 
the wheel moves respectively from 1' to 2' t:hen 3'. 


As the wheel moves 
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from position 1 , a shear surface confining an incremental soil wedge forms 
with an active Rankine zone aibici, a radial shear zone with log-spiral 
center at a*, and a passive Rankine zone. As the wheel continues to move 
(roll), the active Rankine zone grows and new shear surfaces form. Evidence 
for such shear surfaces is shown in Fig. 2-1; shear surfaces have also been 
observed in model tests (Hovland and Mitchell, 1972). The active Rankine 
zone and the corresponding soil wedge continue to grow until a limiting 
(maximum) size of wedge is formed. Beyond that point, new limiting (maximum) 
soil wedges will continue to form as the wheel rolls forward. It has been 
found that the maximum soil wedge (see Fig. 7-1) can be defined approximately 
by 

1. An active Rankine zone, a 2 b 2 C 2 , where distance (a 2 b 2 ) equals one- 
half the soil-wheel contact or 1/2 (a 2 d 2 ), with angle ip - 1.2 <J> 
(Meyerhof, 1955) 

2. A radial shear zone with log-spiral center at a 2 

3. A passive Rankine zone with 02 and 63 as indicated (Sokolovski, 

1969; Karafiath and Nowatzki, 1968) 

As the wheel moves from 2 to 3, the limiting soil wedge will deform and 
move out with relative motion, as shown by the arrows in Fig. 7-1. At this 
point, the active Rankine zone will be more compressed and the center' of 
the log-spiral will be located at the intersection of the 02 line, with 
the vertical drawn through C 33 . At this point, with the center of the log- 
spiral directly above C3, motion along this ultimate shear surface will 
stop, since the angle between the direction of motion 6f the active Rankine 
zone and the normal to the shear surface is cp. For the wheel at position 3, 
the limiting soil wedge, constructed with the center of the log-spiral at 
3, has the same size as that constructed foSr the wheel at 2 with the center 
of the log-spiral at a 2 - 
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To investigate how well the soil wedge selected by the procedure out- 
lined above describes the actual soil wedge, a 24-inch (61 cm) diameter 
spherical wheel was rolled on Yuma sand, and the extent of the shear zones 
in front of and to the sides of the wheel were measured. Fig. 7-2 gives 
the results of this comparison. Predicted distance of forward shear (the 
distance at which the shear surface exits the slope) is plotted in comparison 
with measured distance of forward shear for different wheel loads. The 
measured distance of forward shear is the distance from the crest of the 
track to the furthest shear surface. The measured distance of side shear 
is also shown. It is the width of the shear zone on the side of the track 
measured from the crest of the track. The selected shear surface geometry 
appears to be adequate. 

Mass of soil wedge : 

Due to mathematical complexities in expressing exactly the mass of 
soil bounded by a log-spiral shear surface, and inherent assumptions in 
applying plasticity theory to the wheel-soil interaction problem, an 
approximation for the mass of soil is developed. 

A longitudinal section of the boundary confining the moving soil, 
and an approximation to that boundary confining approximately the same mass 
of soil are shown in Fig. 7-3a. The ratio of d/d Q is plotted in Fig. 7-3b; 
this ratio was determined graphically. The curve in Fig. 7-3b can be 
expressed by 

^ = 1 + 0.002<f> 1,77 (7-2) 

o 

where <j> is the soil friction angle in degrees. 

From Fig. 7-3a, it follows that 



(7-3) 
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w = track width 

O MEASURED SIDE SHEAR 
• MEASURED FORWARD SHEAR 
A PREDICTED FORWARD SHEAR 



RATIO OF MAXIMUM WIDTH OR 
LENGTH OF SHEAR ZONE TO 
AVERAGE TRACK WIDTH 


Fig. 7-2. Wheel load vs. predicted and measured dimension of soil wedge. 




WHEEL LOAD IN POUNDS 
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VERTICAL 

4 _ 



a) Boundary confining approximately the same volume of soil as the 
actual shear surface. 



i 

Fig. 7-3. Geometry for volume of moving soil. 
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where r = radius of wheel, 0^ = wheel-soil contact angle (see Fig. 7-1). 

Fig. 7-2 suggests that the width of the moving soil wedge can be approximated 
by 1.4 w, where w is the track width. Then, the mass of soil can be expressed 
as 


p wd“ 

m = p V =0.7 ~ - 5 — 
s s s tan 03 


where p = soil mass density and V = soil volume, 
s s 


(7-4) 


Acceleration of Moving Soil 

Acceleration of the moving soil wedge is related to its velocity. 
Velocity of the soil wedge is, in turn, related to the velocity of the 
wheel, which is known. It is, therefore, of interest to express the soil 
acceleration in terms of the velocity of the wheel. 

Maximum soil acceleration as a function of average soil velocity : 

Since motion of the soil wedge starts from rest and stops after the 
rolling wheel has passed, the initial and final soil velocities are known 
to be zero. Somewhere in between, the velocity of the shearing soil reaches 
a maximum. Fig. 7-4 shows experimentally determined distance, velocity, 
and acceleration relationships for shearing soil wedges. The plots were 
made from films taken of rolling spheres. Some 200 spheres were rolled on 
Yuma sand carefully prepared to desired densities in large soil cars at 
the U. S. Army Waterways Experiment Station (WES) (Hovland and Mitchell, 
1971), and in some cases it was possible to determine the motion of soil 
wedges. 

The relationships shown in Fig. 7-4 can be approximated theoretically 

by various assumed functions as listed in Table 7-1 and plotted in Fig. 7-5. 

The indicated accelerations, a ; velocities, v : and distances, x ; are 

s s s 

derived from the assumed function by simple ^differentiation and integration. 




Fig. 7-4. Soil movement as a function of time (Hovland and Mitchell, 1971). 



Table 7-1 


Functions Approximating Motion of Soil Wedge 


Assumed Function 


Linear acceleration 


a = — Bit + B2 
s 



Resulting Expressions 


t V a 

s ■ 6 ( 1 - 2 1 r> f 

s s 


v ■ 67(1-7)’ 

s t t a 

s s 


, t ,1 1 t V 

x s * ‘r'l'ir 1 ’. 

s s 


Maximum 

acceleration 



Sinusoidal velocity 


2 . v 

tt r ft i a 

2~ cos [tt (1- — 

s s 


V 

4 . 9 ^ 

s 


v g = Bisin (B 2 +B3t) v 


= x- sin [tt (1 - ~)] V a 
s z t 

s 


X - -5 s - [cos[TT(l - f-)] + 1] V 

s 2 t a 



Step acceleration 


a = constant 
s 


Sinusoidal acceleration a 


a = a sin(Bit+B2) v 
s max s 


t V a 
2 tt sin (2tt— ) — 

s s 


(1 - cos (2 tt— -)) v 
t a 
s 


(t - 2! sin (?7r r _)) v £ 
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In Table 7-1, v = average soil velocity, t = time in which soil moves, and 

3 . S 

Bj , B 2 , and B 3 are integration constants. 

By comparing the experimental curves (Fig 7-4) and the theoretical 
curves (Fig. 7-5) , and the values of maximum acceleration listed in Table 
7-1, the following conclusions can be made: 

1. Since a small difference in the distance vs. time curve (Fig. 7-5) 
can lead to a drastic difference in the acceleration curve, it is 
difficult to interpret the exact acceleration function from an 
experimentally measured distance vs. time curve. 

2. Since the maximum accelerations for. the various functions (Table 

7-1, Fig. 7-5) only vary from 4 v /t to 6.28 v /t , it is not 

3 . s a s 

necessary to know the exact acceleration function, for a first 
approximation . 

3. The curves in Fig. 7-4 show greatest resemblance to the step 

acceleration and sinusoidal velocity functions in Fig. 7-5 for 

which a = 4 v /t and 4.9 v /t , respectively, 
max as as 

In the following development, a = 5 v ft will be used. Deviations 

max a s 

from this value will then be within ± 20 % of the range of functions 
considered. 

Average soil velocity as a function of wheel velocity : 

Fig. 7-1 indicates that in the time the wheel moves from 2 to 3, the 
soil wedge (midpoint of a 2 C 2 ) moves approximately 1/6 that distance in the 
same direction. Therefore, the average soil velocity is 



(7-6) 
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To have a useful expression for a g it is necessary to eliminate t from 
equation (7-6). Observations of wheel-soil interaction show that the soil 
at ai , Fig. 7-1, begins to move before the wheel reaches position 1. With 
the very limited experimental data available on this problem, it will be 
assumed that the soil moves in the time the wheel moves twice the distance 
from position 1 to 3 in Fig. 7-1. Then we can say that t = B^r/v, and 


A v 2 

6 0 T r 


(7-7) 


It is now possible to compare experimental maximum soil accelerations 


shown in Fig. 7-4 and theoretical values determined using equations (7-6) 


and (7-7). This comparison is shown in Table 7-2. 


Table 7-2 

Comparison Between Experimental and Theoretical Accelerations 


Identification 

Soil car Slope Sphere 
° Diam. 



Experimental Theoretical 


V 

(cm/sec) 

t 

s 

(sec) 

a 

s 

(cm/sec 2 ) 

a g (Eq. 7-6) 
(cm/sec 2 ) 

a (Eq.7-7) 
s 

(cm/sec 2 ) 

19 

.25 

67 

63 

51 

27 

.085 

320 

265 

102 

106 

.085 

1060 

1040 

1675 * 


In Table 7-2, the theoretical accelerations computed using equation (7-6) and 


experimental values for t compare well. The theoretical accelerations 

computed using equation (7-7) and the assumed value for t compare less well. 

♦ s 


Soil Inertia Resistance to a Rolling Wheel 


Combining equations (7-4) and (7-7) gives 

0.58p wd 2 v 2 

ft . S / 

I tan 03 B^r 


(7-8) 
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This inertia force can be expressed as force per unit area of wheel to soil 
contact, analogous to unit bearing capacity, by dividing by the contact 
area, 0 T rw, 


0.58 p d 2 2 

-5 ( _Z_) 

tan0 3 ^r' 


(7-9) 


A convenient dimensionless expression results by dividing both sides by p gw 

s 

and substituting for d, equations (7-2) and (7-3) 


p s gw 


(1 + 0.002 <j> 


13.8 tan 03 


1.77 2 2 

)_ v 

wg 


(7-10) 


or 


where 


q I 

p s gw ~ 



(7-11) 


fi = (1 + 0.0020 1 - 77 ) 2 

P 13.8 tan0 3 

As an aid in solving equation (7-11), values of 8 are presented in Fig. 7-6. 
As a approaches <p, 0 3 approaches zero and 6 approaches infinity. Since such 
values of 8 are unreasonable. Fig. 7-6 is considered valid for a < <j> - 5° 
only; i.e., the solid portion of the curves only. 


SCHURING'S ANALYSIS 

It is also possible to study the effects of inertia forces by consider- 
ing the flotation provided a moving wheel by soil inertia and soil static 
forces (bearing capacity) . It is to be noted that while the volume of the 
soil wedge is not. a function of cohesion (previous section), bearing 
capacity is a function of both cohesion and friction. Thus, at this point 
cohesion, c, is introduced. 

Schuring (1968) studied inertia forces by moving plates at various 

t 
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velocities and inclined at various angles to a bed of soft clay. He plotted 
his results as shown in Fig. 7-7. This plot is analogous to the familiar 
Moody diagram in fluid mechanics; the ratio plotted on the vertical axis 
is similar to the drag coefficient and the ratio plotted on the horizontal 
axis corresponds to Reynolds number. Schuring concluded that: 

1. For ordinary vehicular speeds, inertia effects are likely to be 
insignificantly small. 

2. "For special applications, such as dynamic testing and aircraft 
landing, inertia of accelerated soil may be the most important 
factor. " 

Wheel velocities which lead to either insignificant or significant 
(static range and dynamic range, respectively) inertia forces can be 
conveniently separated by plotting data as shown in Fig. 7-7. In the 
static range, where inertia forces are negligible, data plot along a 45° 
line. In the dynamic range, where inertia forces predominate, data 
become asymptotic to the horizontal axis. Thus, this approach places 
bounds on velocities for which inertia forces can be neglected and for 
which they predominate. 

This ingenious approach appears to have two shortcomings, which Schuring 
also recognized. First, although the method is good for separating - 
negligible and predominating inertia forces, it does not provide a basis 
for predicting performance for specific conditions, particularly in the 
transition range. Further, data for a c - 4> soil cannot be accounted for. 

(4> represents an additional parameter, and the frictional force is further 
complicated by being a function of the normal force (see Schuring, 1968)). 

The theoretical analysis presented earlier in this paper attempts to 
overcome the first shortcoming. A concept pf equivalent cohesion pre- 
sented below attempts to overcome the second shortcoming. 



LIFT 
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EQUIVALENT COHESION 

While a general relationship between cohesion, c, and friction, <p, 
is not possible (c and <j> being independent soil parameters) , it is always 
possible to express the total effect of the two by one strength parameter 
for any one type of problem. Soil resistance to a rolling wheel can be 
considered a bearing capacity type of problem. (It has been demonstrated 
that bearing capacity theory can be used as a basis for analyzing soil 
resistance to freely rolling spherical wheels and spheres (see Chapter 5 
and Hovland and Mitchell, 1971). It is believed that phenomenologically 
bearing capacity theory can provide the basis for analyzing soil resistance 
of rolling wheels in general). Thus, for this specific problem, an 
equivalent cohesion is developed below. 

The general bearing capacity equation (5-7), was adapted to the rolling 
sphere problem (Hovland and Mitchell, 1971) 


= 0.188 N + 1.1 


N + 0.55 - 
cs w 


where 


q = soil bearing capacity 

N , N , N = bearing capacity factors (Meyerhof, 1951) 
y c q 

N = 0.25 N 
ys y 

N - 0.50 N 
cs c 

N = N tan<f> + 1 
qs cs 

z = sinkage 

y = soil unit weight = p g 
s s * 

If a soil (x) is to give the same resistance as a soil (1) , 


\ [.188N + 1.1 /— \ N +0.55 {-\ N 1 

wTs/ x I Ys \wy S J cs \w/ qsj x 

) [" . 188N + 1.1 l-~ -) N + 0.55 (-) N 

w Ys Ji L Ys \ W Y S / cs ^Ji 
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Rearranging equation (7-13) and substituting N = 0.25 N , N = 0.5 N , 

jS Y CS C 

and N = N tan<{> + 1 leads to an expression for c /c 
qs cs x ^ 


c 

X 

c 

1 



+ 0.085 13 


N 


cx 


N 


N 


cx 


(7-14) 


In this equation geometrical parameters are held constant ( w x = = w) , 

and Y = Y = Y* Ac-* soil can be treated as an equivalent c soil by 
sx 'si ’ e ■ 

letting <f> x in equation (7-14) equal zero 


e 

wY r 


■[ 


tan<() + 


wyJ N 


N 


ci 


+ 0.085 


cx 


N 

_XL 

N 

cx 


(7-15) 


where N and N are determined for * = d> , and N is determined for * = 0 . 
ci Y'l i cx Y 

Using equivalent cohesion as determined from equation (7-15), experi- 
mental data for any c - <p soil can be analyzed on the type of plot shown in 
Fig. 7-7. Comparing theory and experimental data on this type of plot has 
certain important implications, which are discussed subsequently. 


DISCUSSION 

Theory Compared with Schuring's Data 

To compare predictions of the theoretical analysis presented in this 
paper with Schuring's data, load-soil interaction must be expressed for 
comparable geometries. The geometries of a rolling wheel and Schuring’s 
plate are illustrated in Fig. 7-8. The geometries are comparable, and 
the wheel geometry is defined, if the length and orientation of the chord- 
length for a wheel-soil contact is the same as the length, Z, and the 
orientation of the plate. 

f 



96 


WHEEL 



PLATE 



Fig. 7-8. Comparable geometries for wheel and plate. 


t 
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The axes of Fig. 7-7 can now be expressed in terms of developed theory 
using equation (7-11) and Fig. 7-8: 


lift 

inertia force 


R y * (q + qj) 
p s *^- - -17- 


Bcosot' , and 


inertia force _ v *~s _ ^1 
cohesion force c Be 

e 


where a' = inclination of plate; SL = length of plate in the plane of the 
paper, equal to the chord distance of 0^; and R^= lift, equal to the vertical 
component of the total bearing capacity. In working out the parameters of 
these axes, q^ is computed from equation (7-11), q is computed from equation 
(7-12), assuming a' = slope angle in determining the bearing capacity 
factors (Meyerhof, 1951), and c g is computed from equation (7-15). 

Ordinates and abscissas were computed as described above for the soil 
conditions of Schuring's tests. His experimental curves are compared with 
theoretical predictions, curves A and A' , in Fig. 7-9. While for Schuring's 
experimental curves, the separation in the static range is caused by a 
change in cohesion, c, the separation for the theoretical curves A and A' 
is primarily caused by the change in a' from 15° to 30°. Schuring's data 
for c = 165 psf (8.1 kN/m 2 ) are very close to the theoretical curves. For 
■>the data for c = 35 psf (1.7 kN/m 2 ), the experimental ordinate is more than 
twice the theoretical ordinate for high velocities. Since the mechanisms 
of soil shear under a rolling wheel (theoretical analysis) and a bulldozing 
plate (Schuring's data) are somewhat different, a close agreement for this 
comparison should not be expected. 

The Influence of changes in soil parameters on inertia forces is also 
illustrated in Fig. 7-9. If <j>i is increased from zero to 10°, 20° and 30°, 


*In Schuring's analysis, SL Z = area of plate (verbal communication with Mr. 
Schuring) . 
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keeping all other soil parameters in curve A constant, c £ increases from 
165 psf to 250, 420, and 780 psf (8.1 kN/m 2 to 12.2, 20.6, and 38,2 kN/m 2 ) 
respectively, and we obtain curves B, C, and D. Therefore, in the dynamic 
range, curves A, B, C, and D separate with an increase in <j). If, instead, cx 
is increased, keeping <|>i = 0, curve A remains unchanged. This behavior is 
consistent with plasticity theory, which predicts that the volume of a 
soil wedge is a function of <p but is not a function of c. 

Theory Compared with Data from Rolling Sphere Tests 

An experimental comparison can also be made for spheres rolled on 
Yuma sand (Hovland and Mitchell, 1971). In this case, experimental values 
for the ordinates and abscissas for the type of plot shown in Fig. 7-7 are 
most conveniently determined as R/(p £ 2 v 2 ) and (p v 2 )/c , respectively. 

S S 6 

Dynamic equilibrium of a freely rolling sphere (Hovland and Mitchell, 1971) 
(see also Fig. 4-1 and equation (4-5)) shows that 

R - W[ (sinot - -) + cos 2 a] 1/2 (7-16) 

8 

where R = vertical force or lift, W = sphere weight, and a = sphere 
acceleration. Since W, a, and a are known, R can be calculated from 
equation (7-16). The sphere-soil contact area is % z =0.393 w 2 (Hovland 
and Mitchell, 1971), where w = track width. Theoretical values for the 
ordinates and abscissas are determined as 8(q + q^/q^) and 9j/(^ c e )> 
where q^. is computed from equation (7-11) , q is computed from equation 
(7-12), and c g is computed from equation (7-15). Note that for a freely 
rolling sphere, the forces are vertical as determined, and the cosa 
correction of the ordinate values is not necessary. Table 7-3 lists the 
type of data compared and Fig. 7-10 shows the comparison- 
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Table 7-3 

Data £or Spheres Rolled on Yuma Sand 


Curve 

on 

Fig. 7-10 

Soil 

Cars 

Experimental Data 


Basis 

for Theoretical 

Curves 


Range 

of 

a° 

Range of 
D (cm) 

a ° 
ave 

w 

D 

ave 

(cm) 

p s 

gn-sec 2 
^ cm 4 2 

4>i° 

— 

c 

e 

(gm/cm 2 ) 

A 

1,2 

15-25 

.45-1.00 

20 

.68 

. 00150 

37 

0 

12 

B 

3 

10-35 

.30-. 65 

20 

.40 

.00165 

42 

2.1 

36 

C 

4,8,9 

10-35 

.10-. 47 

20 

.33 

.00150 

37 

10.5 

85 

C 

5,6,7 

10-35 

.32-. 87 

20 

.61 

.00128 

32 

4.8 

28 


If Figs. 7-9 and 7-10 are compared, It will be noted that the data for 
Yuma sand (32° < < 42°) are approximately averaged by Schuring's c = 165 

psf (8.1 kN/m 2 ) line, for which (|> * 0. This suggests that, in the static 
range, data plotted in terms of c g tend to converge to a narrow zone. 

While it has been argued that a change in cohesion has a minor effect, 
curves A and C in Fig. 7-10 suggest that a change in w/D can have a much 
larger effect. However, the curves in Fig. 7-10 provide only a general 
comparison between theory and experimental data. Note that while 
the curves were determined for an average w/D, for the experimental 
data, w/D changed from data point to data point. 

A more direct comparison is shown in Fig. 7-11, where the actual w/D 
ratio for each point was used in computing the theoretical values. In 
general, the comparison is good. The data in both Figs. 7-10 and 7-11 
particularly for cars 1, 2, and 3 show a relatively large amount of scatter. 

Most of the data for the spheres rolled on moist, dense Yuma sand (car 4) 

* 

lie approximately 20% below the 45° line. Most of the data for the spheres 
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rolled on moist, loose Yuma sand (cars 5,6,7) lie approximately 20% above 
the 45° line. 

Since for the soil conditions of these tests (Table 7-3) 3 varies 

approximately from 0.75 to 1.25 (Fig. 7-6) with an average value of 1.00, 

the relative contribution due to soil inertia is approximately the inverse 

of the value of the axes in Fig. 7-11. That is, for R/(p Av z ) = 1, 10, 

s 

and 100, inertia contribution to total soil resistance is approximately 100%, 

10% and 1% respectively. Therefore, to the right of R/(p g Av 2 ) = 10, the 

comparison is indirect evidence for the applicability of bearing capacity 

theory in wheel-soil interaction. For all data to the left of R/ (p Av 2 ) = 

s 

10 in Fig. 7-11, soil inertia forces were significant. This fact is not 
clearly revealed in Fig. 7-10, where the data essentially plot along a 
straight 45° line with no indication of bending into the transition range. 

For low values of the abscissa in Fig. 7-11, the data tend to be above 
the 45° line. Although the data are limited, this suggests that for these 
rolling spheres, the theory underestimated the inertia effects by perhaps 
50%. 

Value of Theory to Experimental Work 

While the agreement between theory and experimental data is generally 
good, although very preliminary, it is equally important that g can be 
determined experimentally. At high velocities (100 mi/hr.), as becomes 
very large, (q^ + q) /q^. approaches one, and the value of the ordinate 
(for example in Fig. 7-10) approaches g. 

It will be recalled that the assumptions used in developing equation 
(7-11) are all incorporated into g. Therefore, as experimental data on g 
become available for various soils, the assumptions can be checked and 
better analysis will be possible. 
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The expression for 3 used in equation (7-11) can be generalized as 


V ip 

3 * 

(0r) 2 w 


(7-17) 


where V = volume of soil, and ip is a function or coefficient relating the 
s a. 

V 2 

acceleration of the moving soil to the velocity of the wheel (a g = ip a g^-) . 

The only unknowns in aquation (7-17) are V and tp . Therefore, experimentally 

S 3 . 

determining 3 is a way of backcalculating V ip . Since ip may be close to 

s a a 

one (see equation (7-7)), experimentally determining 3 is a way of back- 
calculating the volume of moving soil. 

Examples 

To illustrate the practical implications of the theoretical analysis, 
consider first a Cessna 150 attempting a landing on a soft playa or marsh. 

If the gross weight of the airplane is assumed to be 1600 lb. (7300 N) , at 
what velocity will severe and immobilizing sinkage of the wheels begin? (At 
what velocity will the wheel loads exceed the total lift?) Fig. 7-12 shows 
the velocity at which the wheels begin to sink severely into the soil, the 
total lift, and the lift provided by soil inertia forces as a percentage 
of total lift. (Soil and wheel conditions are indicated in Fig. 7-12.) 

The curve shown is computed in accordance with presented theory, and no 
experimental data are available at this time for comparison. 

Consider now the inertial contribution to lift of the Lunar Roving 
Vehicle (LRV) operated on the surface of the moon. Estimates by the 
theoretical analysis are shown in Fig. 7-13. (Soil conditions are as 
indicated.) Inertia forces are shown to be significant at low velocities. 

Note that, on the moon, soil forces which are dependent on gravity are 
reduced by 1/6; therefore, the contributions due to inertia and cohesion 
are relatively six times more important . 1 



lift PER WHEEL, lb 
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CONDITIONS: 

W = 7", i = 5.7", 2 ~ 2" 

4 > -- 10 ° 

C = 500 psf = 24.5 kN/m* 

>S = 100 pcf « 16.0 kN/m 3 

VELOCITY, v , Km/hr 

0 20 40 60 80 100 120 



Fig. 7-12. Relative importance of soil inertia forces for a light 
airplane landing on a soft clayey soil. 


# 
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CONDITIONS : 

W = 9", 1=6", Z = 0.5" 

4> - 40° 

C - 20 psf = 0.98 kN/m z 
y = 100/6 pcf - 16/6 kN/m 3 

o 


VELOCITY , v , Km/hr 



Fig. 7-13. Relative importance of soil inertia forces for the Lunar 
Roving Vehicle operated on the surface of the moon. 


) 


Qi/(q + q x ), % 
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While in the first example (Fig. 7-12) the wheel load is higher than 
the static bearing capacity, the wheel load for the Lunar Roving Vehicle 
(appr. 35 lb = 160 N on the moon) is much less than the static bearing 
capacity (appr. 118 lb = 540 N on the moon). Therefore, not all the soil 
within the soil wedge defined by theory moves. The relative amount of 
soil that moves can, however, be estimated. 

Radial pressure, a, at a wheel-soil contact is proportional to the 
resistance to shear that is mobilized along a shear surface. Resistance 
to shear is proportional to the normal pressure along the shear surface. 

The normal pressure at any point along the shear surface (see Fig. 7-1) can 
be estimated by adding the pressure due to wheel load and the pressure due 
to soil weight. (The pressure due to the wheel load can be estimated from 
theory of pressure beneath a uniform load on an elastic half-space.) 
Assuming that for an increment of deformation, x, of the soil at the 

wheel-soil contact, shear stresses are mobilized along a related portion 

\ 

of the shear surface, we can estimate and plot resulting radial pressure, 
a, vs. displacement, x. Such a plot suggests that the ratio of displace- 
ment, x, to displacement for full mobilization of shear stresses along 
the shear surface, x , can be considered proportional to the ratio of 
a to bearing capacity, q, squared. That is. 


— 0 ( £ f 

X max W 

Therefore, for the LRV wheel-soil interaction, the magnitude of soil 

* 

inertia forces, which are proportional to the amount of moving soil, are 
reduced by approximately (35/118)^ - (1/11). This leads to the dashed 
line in Fig. 7-13. 
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The above example also clearly illustrates that the performance of 

I \ 

a deformable wheel operating below soil bearing capacity is a different 
and a more complex problem than a rigid wheel, and Fig. 7-13 is intended 
as an example and a rough estimate only. Note that while for a rigid 
wheel, equilibrium is achieved by appropriate sinkage into the soil, for 
a deformable wheel, equilibrium is achieved mainly by appropriate flexure 
of the tire. 

Based on the dashed curve in Fig. 7-13, soil inertia forces are 
likely to contribute from 10 to 20 percent of the total lift to the Lunar 
Roving Vehicle. 

CONCLUSIONS 

A theory has been presented from which inertia forces of moving soil, 
influencing the motion of a rolling wheel, can be estimated. The theory 
is compared with previous experimental analysis and experimental data from 
rolling spheres. These comparisons indicate that the theoretical pre- 
dictions are reasonable, and that predictions specifically of the magnitude 
of inertial effects are within 50% of the measured values. These inertial 
effects are computed from equation (7-11) 

q I v 2 

— — = 8 — 

Pggw wg 

where q^ = pressure due to inertia of moving soil (acting in the same 

direction as static bearing capacity, q) , p = mass density of soil, 

s 

g - acceleration of gravity, w = track width, v = velocity of wheel, and 

■* 

3 is a function of soil friction angle, <j>, and slope angle, a. Values 
of 8 are given in Fig. 7-6. 

A Goncept of equivalent cohesion is introduced, which allows data for 
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and c — <j> soil to be plotted and analyzed as shown in figs. 7-7, 7-9, and 
7-10. The main advantage of this type of plot is that 3 can be determined 
experimentally from tests at high wheel velocities (100 m/hr ** 161 km/hr) . 
Since all assumptions necessary in developing equation (7-11) are incor- 
porated in 3, refinements in values of 3 can be made as additional 
experimental data become available. 

The studies presented show that soil inertia forces can be important 
(5 to 50% of total soil reaction) at even moderate wheel velocities. 
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CHAPTER 8. WHEEL-SOIL INTERACTION ANALYSIS 

SUMMARY 

A summary of the relationships derived in the previous chapters 
is given below. These relationships are related to: 

1. Sinkage of a wheel 

2. Pull that can be developed from a wheel 

3. Mobilization of shear stresses at the wheel-soil contact 

4. Soil inertia effects 

The derivations are all based on the assumption that the wheel- 
soil interaction is two-dimensional in the sense that relevant inter- 
action between soil particles and the wheel is all in the plane of 
rolling of the wheel. The presented relationships were based on two 
fundamental observations: 

1. The line of action of the resultant of radial stresses approxi- 
mately bisects the wheel-soil contact angle, 0^, for all values 
of slip, s. 

2. For developing general wheel-soil interaction theory, all 
parameters operating at the wheel-soil interface are functions 
of slip and position. This led to concepts such as 1) slip at 
a point, 2) the shear stress surface, and 3) the performance 
surface (to be presented) . The general problem requires that 
all functional relationships be visualized in some three- 
dimensional way. 

Summary of Basic Relationships 

For a rigid cylindrical wheel operating in any soil, sinkage can be 
related to the towed force by 




Ill 


N 

/ - tan [0.4 cos' i (l-2 §)) 

y 


(5-4) 


where z - sinkage, D = wheel diameter, N x = motion resistance, and 

Ny = vertical component of N associated with motion resistance (Fig. 4-10). 

Sinkage can be estimated, using bearing capacity theory, from equation 
(5-8) or equation (5-9) by noting that the total bearing capacity, Q, must 
equal, N, giving 

q = f " -f- + s cq [(c+q' tan<j>) N £ + q'] (5-8) 

where q = unit bearing capacity 

q' = z Y /2 or surcharge 
s 

N = resultant of radial stresses 

A = wheel-soil contact area = © T rw 

w = wheel width 

y = soil unit weight 
s 

c = soil cohesion 



Sinkage can also be estimated from 
k 

q = [— + k , ] z n (5-5) 

w 9 

using procedures outlined by Bekker (1969) , and taking into account 
experience gained from the us'te of these procedures. 


i 
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Pull, P, negative or positive, can be determined from a force polygon 
(Fig. 4-10) or from the following equations 

For a soil with both cohesion and friction, where the wheel-soil cohesion 
c >0, and the wheel-soil friction angle 6 > 0, 

£L 


(s) 2 ■ ( cos6 n - 1 [ i+ (=7r) ] 


2 (cose,, -t sine„\ 

\ N W N j ( r / r )Z H Vwr/r/ 


For a pure sand (c « 0) , 

3 


P ^ as 


r 

f a0 del 2 

, r f /r £ J 

f cos 6 


(I) 2 ■ ( cos8 n - 5 Sln6 s) 2 + (Itt) 2 - 


For a pure clay (6 *= 0), 
(I) ■ ( COS0 H 


cos0„ - 77 sin8„ 
N W N 


/Ac \ 2 


For a towed wheel in any soil, 


(I)" • ( COS0 N ' I 8ln0 n) ' 
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The towed force can also be analyzed from 



(4-31) 


Soil inertia forces influencing the wheel can be evaluated from 


q i = g vl 

p Q gw gw 


(7-11) 


where 3 can be estimated from Fig. 7-6. 

In the above equations, 6 and c are functions of position, 0, 

c l 

and slip, s, but the subscripts were omitted to shorten the expressions. 

(The use of 6 and c and the mobilization of shear stresses are discussed 

a 

in detail in Chapter 6.) The ratios "r/r and r^/r^, which turn out to be 
equal, can be determined from Fig. 3-7. The terms in the above 
equations are defined as follows : 


c = wheel-soil cohesion 
a 

g = acceleration of gravity 
P = pull (negative or positve) for any slip 
q^- = soil inertia pressure 

r = radius of wheel 

— = radius to line of action of F (vectorial) 

r^ = moment arm to dR (see Fig. 3-1) 
r^ = moment arm to R (see Fig. 3-1) 
v = translational velocity of wheel 
W = wheel load 

3 “ abbreviation used in solving equation (7-11), see Fig. 7-6 

6 = wheel-soil friction angle 

0 = angle identifying any positions at the wheel soil contact 

(see Fig. 3-1) 
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0 N = angle to the line of action of the resultant of radial (normal) 
stresses 

Cfg = radial stress at the wheel-soil contact 

p = soil mass density 
s 

ANALYSIS PROCEDURE 

Suppose that we are to determine performance of a wheel in a given 
soil at a particular state of slip of the wheel. In solving all forces 
of the force polygon (Fig. 4-10), two separate procedures may be considered: 
1) sinkage known, 2) sinkage unknown. These cases will be described in 
detail below. The written outline follows the block diagram and 
graphical solutions in Figs. 8-1 and 8-2. 

Sinkage Known (Fig. 8-1) 

Given: Wheel load 

Procedure: 

1. Determine wheel-soil strength parameters, 6 and c a , as a 
function of displacement from 3 ring-shear tests, as shown in 
Fig. 6-5. 

2. Select the appropriate strength parameters, 6 and c_, for desired 

a 

slip, using Fig. 6-4 and the test data. 

3. Compute 0^ from equation (5-3), and draw to scale W and 0^. 

4. Assume a value of N and compute F using equations (4-18) and 
(4-19). Draw F perpendicular to the 0^ line, and draw N 
parallel to the 0^ line, from F. 

5. Compare F and N to the required closure. The resultant of F and N 
must reach the line of action. at P, which is a horizontal line drawn 
from the top of W. If the force polygon does not close, a wrong N 

t 


value was used. 
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Fig. 8-1. Solution of the wheel-soil interaction 
force polygon for known sinkage. 
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6. Select a new N and compute F; continue steps 4, 5 and 6 
until the force polygon closes. 

7. Closure of the force polygon determines the pull, P, and all 
other forces associated with the force polygon. Once all 
these forces are known, the required torque, or torque 
associated with F, can also be determined. 

Weakest Link: 

The weakest link or step in the above procedure is probably 

that of relating soil parameters to slip (Fig. 6-4). This requires 

further research. 

Sinkage Unknown (Fig. 8-2) 

Given: Wheel load 

Procedure: 

1. Determine soil strength parameters, (j> and c, and/or coefficients 
associated with equations (5-5) or (5-6) . 

2. Compute q vs. ^ from equations (5-8) or (5-9) or equations 
(5-5) or (5-6) , using also equation (5-3) . 

3. Assume a value of N giving a value for 0^. Draw to scale W and 0^. 

4. Determine wheel-soil strength parameters, S and c , as a function 
of displacement from 3 ring-shear tests, as shown in Fig. 6-5. 

5. Select the appropriate strength parameters, 6 and c , for 
desired slip, using Fig. 6-4 and the test data. 

6. Compute F using equation (4-18) and (4-19) , 'and draw F perpendi- 
cular to the 0^ line Draw N parallel to the 0 N line, from F. 

7. Compare F and N with the required closure. 
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Fig. 8-2. Solution of the wheel-soil interaction 
force polygon for unknown sirikage. 


} 
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8. If the force polygon does not close, select a new N, compute 
a new F, and continue steps 3, 6, 7, and 8 until appropriate 
closure is achieved. 

9. Closure of the force polygon determines all the forces 
consistent with the computed sinkage and wheel-soil contact 
angles. Again, the torque, T, associated with F can now be 
determined . 

Weakest Link: 

The weakest link or step in the above procedure (sinkage 
unknown) is probably that of estimating sinkage (step 2) , either 
using bearing capacity theory or equations (5-5) or (5-6) . This 
requires further research. In addition, step 5, which relates soil 
parameters to slip (Fig. 6-4), requires further research. 

COMPARISONS 

It has not been possible to compare predictions based on the 
procedures outlined above directly with experimental data, because 
the type of displacement dependent wheel-soil strength parameters 
(Fig. 6-5) required have not been found in the literature. It is 
possible, however, to make some general comparisons. 

A common expression for developed pull (Bekker, 1969) is 

P = H - (8-1) 

H = Ac + Wtantj) 

kwD 

where U = traction developed by the wheel, “ motion resistance, k «* a 

sinkage coefficient, and the exponent x varies between 1 and 2. 

! 
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Consider the simple case of a pure sand with 6 = 45° , k = k^ = 3, 
and a wheel with w = 5 and D = 20. Equation (8-1) then reduces to 


P = W - 


w" 

300 


( 8 - 2 ) 


Based on theory presented in this report, the same simple case can 
be studied using equation (4-21) . For the same soil conditions used 
above, and assuming r/r = 1, equation (4-21) reduces to 

P 2 = 2(Wcos0. 7 - Psin8„) 2 - W 2 (8-3) 

N N 

While equation (8-2) is a function of the width of the wheel, 
equation (8-3) is just an equilibrium solution of forces in the plane of 
rolling. Equation (8-3) incorporates also the case of an undeformable 
surface, such as may be the case with a rail. Note that in this case 
0^=0, and equation (8-3) reduces to P/W = 1. This is obviously correct 
for 6 = 45° . 

For a deformable surface, such as a soil, 0„ is a function of W. 

Assume for the sake of this comparison that 0^ increases linearly from 
0 to 20 degrees as W increases from 0 to 300. On this basis, equations 
(8-2) and (8-3) are compared in Fig. 8-3. As shown, the two theories lead 
to similar P vs. W relationships for this simple case. In the use of 
equation (8-2), the trick is to select the proper x. In the use of equation 
(8-3), the trick is to select the proper 0^. 

The influence of the uncertainty in selecting 0^ was also evaluated 
by solving equation (8-3) for 0^ ± one standard deviation, based on data in 
Figs. 4-5 to 4-9. The range suggested by these solutions is shown by the 
cross-hatched area in Fig. 8-3. The same range can also be defined by 
x = 1.87 and x = 1.91. While it may be difficult to even distinguish 

between 1.87 and 1.91 by curve-fitting techniques, this range represents 

! 




approximately 68% of the data shown in Figs. 4-5 to 4-9. The potential 
value of presented theory, rather than theory based on an exponential 
equation, is therefore demonstrated. Solutions can be more accurately 
bracketed by 0„ rather than x. 

Consideration of pull as a function of both slip, s, and position, 0, 
suggest a performance surface as illustrated in Fig. 8-4. Reanalysis of 
data by Leflaive, 1967, demonstrates experimentally the same type of surface, 
as shown in Fig. 8-5. There is a separate P vs. s relationship for each 
value of W. Therefore, P/W by itself, although very frequently given in 
the literature, may not be a very meaningful dimensionless parameter. 

It is also possible to predict a performance surface from theory 
presented in this report. Such a surface (Fig. 8-6) was computed on the 
basis of Leflaive's wheel and sinkage data, and Figs. 6-4 and 6-5. It is 
to be noted that Figs. 8-5 and 8-6 are not comparable; a comparable 
theoretical prediction could only be made using soil strength parameter 
relationships based on actual testing instead of Fig. 6-5. Fig. 8-6 
illustrates, however, the general and realistic nature of the prediction. 

Also, the performance surface of Fig. 8-6 suggests a more plastic soil than 
that of Fig. 8-5, where there is essentially no increase in pull, P, beyond 
s = + . 15. 

CONCLUSIONS AND RECOMMENDATIONS 
Conclusions 

Approximate and general solutions have been developed for a rigid cylindrical 
wheel operating in soil. A detailed summary of these developments was pre- 
sented earlier in this chapter. The solutions deal with 

1. Sinkage of a wheel 

2. Pull that can be developed from a wheel 
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Fig. 8-5. Performance surface, rigid wheel data, 
Leflaive (1967) 
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3. Mobilization of shear stresses at the wheel-soil contact 

4. Soil inertia effects 

The solutions are primarily the result of simple considerations of statics 
and dynamics. The error associated with the solutions can be evaluated, 
and theoretical predictions are likely to be within 15% of the correct value. 

For each of the four phenomena listed above, the developments suggest 
how prediction of wheel-soil interaction can be improved, and these items 
are discussed below. 

Rec ommendat ions 
Sinkage of a wheel ; 

Sinkage, as computed from equation (5-7), is subject to typical 
limitations of bearing capacity theory. Because the wheel-soil interaction 
problem requires an exact solution as opposed to an upper or lower bound, 
if bearing capacity theory is to be used, it must be tailored specifically 
to fit the problem. This is possible, as was demonstrated by modifying 
the bearing capacity equation to the case of rigid spherical wheels (Fig. 

5-3, 5-4, and 5-5). The intermediate cases (3) and (4) in Fig. 5-6a, 
require further research. 

Tests should be conducted where soil conditions are controlled, and 

careful measurements are made of wheel-soil contact angles, 9^, 

and 0„ as these are influenced by the state of slip of the wheel. Then 
N 

the measured wheel load and pull or towed force can be compared to predic- 
tions by the bearing capacity equation, and more realistic shape factors 
and correction factors to the bearing capacity factors. can be developed. 

The observation that the line of action of the resultant of radial 
stresses approximately bisects the wheel-soil contact angle for all values 
of slip, which is substantiated by data in Figs. 4-5 and 4-6, clearly 

I 
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suggest that the mobilization of shear stresses along the wheel-soil 
contact, and the mobilization of radial stresses, are associated with two 
separate mechanisms. It is believed that the resultant of radial stresses, 

N, should be related to bearing capacity for all values of slip. 

Pull that can be developed from a wheel ; 

Further, while much data has been accumulated on performance para- 
meters, much less detailed data is available on wheel-soil contact angles. 
Since it has been demonstrated that, knowing 0^, the wheel-soil interaction 
problem becomes determinate, careful data on 0^ should be collected. With 
the data presented in Figs. 4-5 through 4-9, it is probable that about 68% 
of predictions of pull will be within approximately 15% of the correct 
value. Improved data, adequate to separate out variables associated 
with the states of the wheel and the soil, should reduce the error to an 
acceptable value in probably all cases. 

Mobilization of shear stresses at the wheel-soil contact : 

Further theoretical contemplation should be first undertaken to relate 
soil test displacements to slip at a point, Sq , as was attempted in Fig. 6-4. 
This could lead to a general and rigorous relationship. Contemporaneously, 
shear stress surfaces (Fig. 6-5) should be plotted from experimental data 
whenever possible to further our understanding of the general problem.* 

Soil inertia effects : 

The developments in Chapter 7, and particularly Figs. 7-9 and 7-10, show 
that 8 in equation (7-11) 

q t - v 2 

P s g» wg 

can be experimentally determined. Since all assumptions necessary in de- 
veloping equation (7-11) are incorporated in 8, our dependence on these 
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assumptions need only be temporary. Tests conducted by running a wheel 
at relatively high velocity on a soil, and measuring the lift provided 
the wheel by the soil, will allow determination of 0 for desired soil 
conditions from the type of plots shown in Figs. 7-9 and 7-10. 
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wheel-soil contact area 



linear acceleration of wheel 
linear acceleration of soil 
maximum acceleration of soil 
integration constants 

width of footing in the bearing capacity equation 
soil cohesion 

wheel-to-soil cohesion or adhesion 

equivalent cohesion 

cohesion of soil (1) 

cohesion of soil (x) 

wheel diameter 

distance revolved 

distance traveled 

dimension defining triangle base proportional to d Q 
dimension of cone (active Rankine zone) underneath wheel 
base of natural logarithm 

vectorial sum of shear stresses at wheel-soil contact 
algebraic sum of shear stresses at wheel-soil contact 
x component of F 
y component of F 
acceleration of gravity 

traction developed by wheel, equation (8-1) 
mass moment of inertia of wheel 



slip = 1 - DT/DR 

cohesion parameter in equation (5-5a) 

friction parameter in equation (5-5a) 

length of wheel-soil contact in direction of mptiopi 

mass of wheel 

mass of soil set in motion by wheel 
Meyerhof (1951) bearing capacity factors 
bearing capacity factors for rolling spheres 
vectorial sum of normal stresses at wheel-soil contact 
x component of N (motion resistance) 
y component of N 

pull developed by the wheel (negative or positive) 

towed force 

total bearing capacity 

unit bearing capacity 

unit resistance due to soil inertia 

surcharge pressure 

vectorial sum of all stresses at the wheel-soil contact 

algebraic sum of all stresses at the wheel-soil contact 

motion resistance, equation (8-1) 

soil inertia force 

x component of R 

y component of R 

wheel radius 

wheel effective radius 

moment arm to any dR, but r^ ^ f (0) 

moment arm to any dR from wheel center 

moment arm to R from wheel center 
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r radius of shear ring test, chapter $ 

r w radius of wheel, chapter 6 

r wheel radius as a function of z 

z 

r radius, or distance from wheel center, to actiop of P 

s slip = (DR - DT) / (DR + DT) 

Sg slip, s, as a function of 0 

s s s shape factors used in bearing capacity theory. 

Y c q 

T input torque 

t time 

t time of soil movement 

s 

t time increment for ring shear test 

t time increment for wheel 

w 

u angular acceleration of wheel 

V volume of moving soil 

s 

v linear velocity of wheel 

v average soil velocity 

v g velocity of soil 

v , velocity of wheel center w.r.t. x 

C/ X 

v^ c velocity of point p w.r.t. c 

v p ^ x velocity of p w.r.t. coordinate system, x, fixed ff» ppfl 

W wheel load or weight 

w track width 

x, y, z coordinates 

x g distance of soil displacement 

x w distance of wheel displacement 

y vertical distance to Instantaneous center of rotation 

z s inkage 


/ 



slope angle 
inclination of plate 

abbreviation used in solving equation (7-1) 

angle defining Meyerhof's free surface 

wheel-soil friction angle 

5 as a function of 6 and s 

unit weight of soil 

unit weight of sphere or rock 

angle identifying any position at the wheel-soil contact 

angle defining the line of action of F 
total wheel-soil contact angle 

angle to the line of action of the resultant of normal (radial) 
stresses 

contact angle associated with ring sheer test, chapter 6 

contact angle associated with wheel, chapter 6 

wheel-soil contact angle forward of vertical centerline 

wheel-soil contact angle rear of vertical centerline 

angle of intersection of slope and upper epd of passiye Rankine 
zone 

angle of intersection of slope and shear surface 
angle defining the direction of any dR in a c - $ soli 
soil friction angle, or angle of shearing resistance 
friction angle of soil (1) 
friction angle of soil (x) 

angle defining frictional resistance at wheel-soil contact 

function or coefficient relating soil acceleration tp wheel 
velocity 
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p g soil mass density 

a radial pressure at wheel-soil contact 

a Qz as a function of 6 and direction perpendicular to the paper, a 

x shear stress at the wheel-soil contact 

t q shear stress as a function of 9 and z 

t)z 

e soil strain associated with wheel 

w 

e soil strain associated with a ring shear test 

At soil displacement associated with a ring shear te?t 

Aw soil displacement associated with wheel 

Az^ depth of soil strained in the ring shear test 

Az depth of soil strained by wheel 

w 

oj angular velocity of wheel 



